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Dust plays an increasingly important role in astrophysics. Historically, 
dust was first recognized for its obscuring effects, and the need to correct 
observed intensities for attenuation by dust continues today. But with the 
increasing sensitivity of IR, FIR, and submm telescopes, dust is increasingly 
important as a diagnostic, with its emission spectrum providing an indica- 
tor of physical conditions, and its radiated power bearing witness to star 
populations of which we might otherwise be unaware. Finally, and most fun- 
damentally, dust is now understood to play many critical roles in galactic 
evolution. By sequestering selected elements in the solid grains, and by cat- 
alyzing formation of the H2 molecule, dust grains are central to the chemistry 
of interstellar gas. Photo electrons from dust grains can dominate the heating 
of gas in regions where ultraviolet starlight is present, and in dense regions the 
infrared emission from dust can be an important cooling mechanism. Finally, 
dust grains can be important in interstellar gas dynamics, communicating 
radiation pressure from starlight to the gas, and providing coupling of the 
magnetic field to the gas in regions of low fractional ionization. 

We would like to understand these effects of dust in the Milky Way, in 
other galaxies, and as a function of cosmic time. These lectures are organized 
around topics in the astrophysics of dust in the the Milky Way, as this is our 
best guide to understanding and modelling dust long ago and far away. 

1 Introduction to Interstellar Dust 

We begin with a brief review of some of the observational evidence which 
informs our study of interstellar dust. 

1.1 Interstellar Extinction 

Through study of open star clusters in the Galaxy, Trumpler (1930) found 
that distant stars were dimmed by something in addition to the inverse square 
law, and concluded that interstellar space in the galactic plane contained "fine 
cosmic dust particles of various sizes ... producing the observed selective ab- 
sorption" . Over the past 7 decades we have built on Trumpler's pioneering 
study, but many aspects of interstellar dust - including its chemical compo- 
sition! - remain uncertain. Let us therefore begin by reviewing the different 
ways in which nature permits us to study interstellar dust. 
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Fig. 1. Extinction at wavelengtli A, relative to the extinction at / = 900 nm, 
as a function of inverse wavelength A~^, for Milky Way regions characterized by 
different values of Ry = Ay / E{B~V), where j4s is the extinction at _B = 4400 ^m, 
Av is the extinction at V = 5500 /^m, and the "reddening" E{B — V) = Ab ~ 
Ay- Note the rapid rise in extinction in the vacuum ultraviolet (A ^ 0.2 /im) for 
regions with Rv ^ 4. The normalization per H nucleon is approximately Aj/Nu ~ 
2.6 X 10~ cm /H. The silicate absorption feature f ill. 5^ at 9.7 ^m and the diffuse 
interstellar bands ( til. 61 are barely visible. 

Trumpler analyzed the interaction of light with interstellar dust, and this 
remains our most direct way to study interstellar dust. We use stars as "stan- 
dard candles", and study the "selective extinction" - or "reddening" - of 
starlight by the dust. With the assumption that the extinction (= absorption 
+ scattering) goes to zero at wavelengths A — > oo, and observing the star at 
sufficiently long wavelength (to, in effect, determine its distance) one can de- 
termine the attenuation of the starlight by dust as a function of wavelength. 
Because atomic hydrogen absorbs strongly for hv > 13.6 eV, it it is only pos- 
sible to measure the contribution of dust to the extinction at hv < 13.6 eV, 
or A > 912A. A typical "extinction curve" - the extinction as a function 
of wavelength or frequency — is shown in Figure ^ showing the rapid rise 
in extinction in the vacuum ultraviolet. Observed extinction curves vary in 
shape from one linc-of-sight to another, but appear to approximately form a 
one-parameter family (Cardelli et al. 1989); the parameter is often taken to 
be the ratio Ry = {Ab — Ay) /Ay, where Ab and Ay are the extinctions 
measured in the B (4400A) and V (5500A) spectral bands. A parametrization 
of the extinction curve was provided by Cardelli et al.; the curves in Fig.^ 
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were calculated using a more recent prescription by Fitzpatrick (1999), with 

the extinction in the infrared following Drainc (1989b). 

We will discuss dust grain optics below, but it is clear that if the dust 
grains were large compared to the wavelength, we would be in the "geometric 
optics" limit and the extinction cross-section would be independent of wave- 
length. Therefore the tendency for the extinction to rise even at the shortest 
wavelengths where we can measure it tells us that grains smaller than the 
wavelength must be making an appreciable contribution to the extinction at 
all of observed wavelengths. As we will see below, "small" means (approxi- 
mately) that 27ra|m — 1|/A ^ 1, where m(A) is the complex refractive index. 
Thus if |m — 1| « 1 at A = 0.1 /im, 

• we must have large numbers of grains with a ^ .015 /xm. 
1.2 Scattering of Starlight by Dust Grains 

When an interstellar cloud happens to be unusually near one or more bright 
stars, we have a "reflection nebula" , where we see the starlight photons which 
have been scattered by the dust in the cloud. The spectrum of the light 
coming from the cloud surface shows the stellar absorption lines, showing that 
scattering rather than some emission process is responsible. By comparing 
the observed scattered intensity with the estimated intensity of the starlight 
incident on the cloud, it is possible to infer the albedo of the dust - the ratio 
of scattering cross section to extinction cross section. The result is that in 
the optical the interstellar dust mixture has an albedo a; « 0.5 scattering 
is about as important as absorption - and the grains are somewhat forward 
scattering, with (cos 6) « 0.5. Rayleigh scattering by particles small compared 
to the wavelength has (cos^) « 0, so this tells us that 

• the particles dominating the scattering at A w 0.6 jum 
have a ^ X/2tt 0.1 /xm. 



1.3 Polarization of Starlight 

The polarization of starlight was discovered in 1949 (Hall 1949; Hall & Mike- 
sell 1949; Hiltner 1949a,b). When it was realized that the degree of polariza- 
tion tended to be larger for stars with greater reddening, and that stars in a 
given region of the sky tended to have similar polarization directions, it was 
obvious that the polarization is produced by the interstellar medium: light 
propagating through the interstellar medium becomes linearly polarized as a 
result of preferential extinction of one linear polarization mode relative to the 
other. The polarization percentage typically peaks near the V band (5500A), 
and can be empirically described by the "Serkowski law" (Serkowski 1973): 

p(A) « p(Ai„ax) exp[-i4:in^(A/Amax)] , (1) 
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Fig. 2. Linear polarization in the ultraviolet measured by Clayton et al. (1992). 
The solid line is the "Serkowski law" fit to the data. Figure from Clayton et al 
(1992). 

with Aniax ~ 5500A and K « 1.15. This "hnear dichroism" of the interstellar 
medium is due to dust grains which are partially aligned by the interstellar 
magnetic field. The peak polarization Pmax is found to fall within an envelope 
< Pmax ^ 0.09(i!/(i? — l^)/mag), or < pv ~ O.OStv; the maximum values 
are presumed to arise on sightlines where the magnetic field is uniform and 
perpendicular to the line-of-sight. While the "Serkowski law" was put forward 
as a fit to the observed polarization at 0.3 /im ^ A ^ 1 /xm, it turns out to give 
a surprisingly good approximation to the measured linear polarization in the 
vacuum ultraviolet (Clayton et al 1992, Wolff et al 1997) as seen in Figure|21 
The mechanism responsible for the grain alignment remains a fascinat- 
ing puzzle, which we will discuss below. Independent of the grain alignment 
mechanism, however, we can infer the sizes of the interstellar grains respon- 
sible for this polarization by noting that the extinction rises rapidly into the 
UV whereas the polarization drops (Kim & Martin 1995). This can be un- 
derstood if the grains responsible for the polarization have sizes a such that 
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Fig. 3. Infrared extinction, relative to extinction at / = 900 nm, showing the strong 
9.7 /im and 18 fim silicate features, the 3.4 /xm aliphatic C-H stretch, and a weak 
unidentified DIB at 1.32 ^tm (Joblin et al. 1990). 

a ~ (Amax/27r) « 0.1 /im: then as one proceeds into the UV one moves toward 
the "geometric optics" limit where both polarization modes suffer the same 
extinction, so the polarization goes to zero: 

• The extinction at A w 0.55 fim has an appreciable contribution from 
grains with sizes a w 0.1 /im which are nonspherical and substantially 
aligned. 

• The grains with a ^ 0.05 /im which dominate the extinction at A ^ 0.3 /xm 
are either spherical (which seems unlikely) or minimally aligned. 

1.4 Spectroscopy of Dust: The 2175A Feature 

Of what is interstellar dust composed? One may look for spectroscopic clues in 
the extinction. The extinction curves in Fig. ^ show a conspicuous extinction 
feature at A^^ = 4.6 /im^^, or A = 2175A. The feature is well-described by 
a Drudc profile. The central wavelength is nearly identical on all sightlines, 
but the width varies significantly from one region to another (Fitzpatrick & 
Massa 1986). 

The strength of this feature implies that the responsible material must be 
abundant (Draine 1989a): it must be made from H, C, N, O, Mg, Si, S, or 
Fe. Small graphite grains would have a strong absorption peak at about this 
frequency (Stecher & Donn 1965; Draine 1989a), due to tt ^ tt* electronic 
excitations in the sp^-bonded carbon sheets. Since the carbon skeleton of 
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polycyclic aromatic hydrocarbon (PAH) molecules resembles a portion of a 
graphite sheet, such molecules also tend to have strong electronic transitions 
at about this frequency. It therefore seems likely that the 2 175 A feature is 
due to some form of s^)^ -bonded carbon material. 



1.5 Spectroscopy of Dust: The Silicate Features 

There is a conspicuous infrared absorption feature at 9.7 fim, shown in Figure 
121 Silicate minerals generally have strong absorption resonances due to the 
Si-0 stretching mode near 10 /im, and it seems virtually certain that the 
interstellar 9.7 /im feature is due to silicates. This conclusion is strengthened 
by the fact that the 10 /im emission feature is seen in the outflows from 
oxygen-rich stars (which would be expected to condense silicate dust) but 
not in the outflows from carbon-rich stars. The interstellar 9.7 fim feature 
is seen both in emission (e.g., in the Trapezium region in Orion [Gillett, 
Forrest, et al. 1975]) or in extinction in the interstellar medium (Roche & 
Aitken 1984). 

Crystalline silicate minerals generally have sharp features in their 10 /im 
absorption which are not seen in the broad interstellar 10 /im feature, leading 
to the conclusion that interstellar silicates are probably amorphous. 

Near 18 /im warm interstellar dust shows another emission feature, which 
is attributable to the Si-O-Si bending mode in amorphous silicates. 




Fig. 4. Extinction at wavelength A, relative to the extinction at / = 900 nm, show- 
ing some of the diffuse interstellar bands. 
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1.6 Spectroscopy of Dust: Diffuse Interstellar Bands 

The 3 features at 0.22 /zm, 9.7 /im, and 18 /im are by far the strongest features 
seen in diffuse interstellar dust. There are, in addition, numerous weaker fea- 
tures in the optical known as the "diffuse interstellar bands" , or DIBs. Figure 
0]shows the extinction for 1.5 /zm^^ < A^^ < 1.75 fim~^, with several conspic- 
uous DIBs present, most notably the DIB at 0.5780 /im. The strongest DIB 
falls at 443. Onm. Jenniskens & Desert (1994) report a total of 154 "certain" 
DIBs in the interval 0.38 — 0.868 /im, plus another 52 "probable" detections. 
DIBs were discovered 80 years ago (Heger 1922) and their interstellar nature 
was established 68 years ago (Merrill 1934). 

It is embarassing that Nature has provided astrophysicists with this wealth 
of spectroscopic clues, yet as of this writing not a single one of the DIBs has 
been convincingly identified! It seems likely that some of the DIBs may be 
due to free-flying large molecules; this hypothetis has received support from 
high resolution spectra of the 5797A feature (see Figure showing intrinsic 
ultrafine structure (Sarre et al. 1995; Kerr et al. 1998). 

1.7 Spectroscopy of Dust: The 3.4/xm Feature 

There is a broad absorption feature at 3.4 /im which is almost certainly due 
to the C-H stretching mode in hydrocarbons. A recent study by Pendleton 
& AUamandola (2002) concluded that hydrocarbons with a mixed aromatic 
(ring) and aliphatic (chain) character provided a good fit to the observed 
interstellar absorption, including the 3.35 — 3.53 /im region. This included 
hydrocarbon films deposited following laser-ablation of amorphous carbon in 
Ar, followed by exposure to atomic H (Mennella et al. 1999) or from a weakly- 
ionized plasma produced by laser-ablation of graphite in hydrogen (Scott & 
Duley 1996; Duley et al. 1998). 

1.8 Spectroscopy of Dust: Ice Features 

In dark molecular clouds a number of additional absorption features are seen, 
most notably a strong band at 3.1 /im which is attributed to the O-H stretch- 
ing mode in H2O ice. However, the 3.1 /im feature is not seen on sightlines 
which are confined to diffuse interstellar clouds (Gillett, Jones et al. 1975). 

When a strong 3.1 /tm feature appears in absorption, a number of other 
absorption features are also seen, including features due to CO (4.67 /im), 
CH3OII (3.53 /tm), and CO2 (15.2 /im). The shape of the 3.1 /im II2O feature 
is indicative of the type of ice and the impurities present in it. The relative 
strengths of the various features indicate that II2O is the dominant "ice" 
species, with NII3, CO, CH3OII, and CO2 as secondary constituents. 
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Fig. 5. Fine structure of tiie 5797A DIB toward ( Opli, ( Per, and /i Sgr (from 
Kerr et al. 1998). Note in particular tlie very narrow feature labelled 0. 

1.9 Spectroscopy of Dust: PAH Emission Features 

A wide variety of galactic objects, including planetary nebulae, HII regions, 
photodissociation fronts, and reflection nebulae, have strong infrared emis- 
sion in the 3-13 //m region. Most of the radiated power emerges in 5 broad 
infrared bands, at 3.3, 6.2, 7.7, 8.6, and 11.3 /im. In Figure IHl we show the 
4-15 /im emission observed from the reflection nebula NGC 7023 (Cesarsky 
et al. 1996). These features are seen in many other objects as well (see, e.g., 
Boulanger et al. 1998). The emission is quite strong: a significant fraction of 
the starlight energy incident on the reflection nebula is reradiated in these 
infrared emission bands - the fraction can be as large as 10-20% depending 
on the spectral type of the star ~ so the particles responsible for this emission 
must be quite abundant, as we will discuss further below. 

The emission is thought to result from vibrational modes of polycyclic 
aromatic hydrocarbon (PAH) molecules. PAHs can be thought of as a section 
of a graphite "sheet" with the C atoms in a 2-dimensional hexagonal lattice, 
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Fig. 6. Emission spectrum of the reflection nebula NGC 7023 measured by ISO- 
CAM (Cesarsky et al. 1996). 



with the H atoms attached to C atoms at the edges of the lattice. Three 
examples of PAHs are shown in Figure [T] When H atoms are attached to the 
edge of an aromatic ring, 

• the C-H stretching mode tends to be at 3.3 /xm, 

• the C-H in-plane bending mode is near 8.6 //m 

• the C-H out-of-plane bending mode tends to fall at 11.3, 11.9, 12.7, or 
13.6 /im, depending on whether there are one, two, three, or four adjacent 
C-H bonds ("solo-", "duo-", "trio-", or "quartet-" C-H). The observed 
spectra suggest approximately equal fractions of solo, duo, or trio C- 
H, and little quartet C-H, consistent with what one would expect for a 
mixture of large symmetric compact PAHs (Stein & Brown 1991). 

• The carbon skeleton has C-C-C bending modes near 6.2 and 7.7 /xm 

The precise position of the bands, and their relative strengths, vary from 
one individual PAH to another, and upon their state of ionization, but the 
observed spectra are in good agreement with what one might obtain from 
a mixture of neutral and ionized PAHs (AUamandola, Hudgins, & Sandford 
1999). 
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Fig. 7. Three examples of PAH molecules. 
1.10 IR and FIR Emission 

The energy which grains absorb from starhght is reradiated at longer wave- 
lengths, mainly into the far-infrared. We can estimate the emission spec- 
trum of the Galaxy using data from all-sky surveys by the IRAS and COBE 
satellites, and pointed observations by the Japanese IRTS satellite. At high 
galactic latitudes we know that most of the dust is illuminated by the local 
average starlight background. Theoretical studies of dust grain heating by 
starlight suggested that a ^ 0.01 /J,m interstellar grains would be heated to 
temperatures 15K ;$ T 20K by the diffuse starlight (see, e.g., Draine & 
Lee 1984). 

Because there is little dust at high galactic latitudes, the infrared sur- 
face brightness is low, but from the correlation of infrared surface brightness 
with 21 cm emission one can extract the infrared emission per H atom. This 
has been done using photometry in the 4 IRAS photometric bands, the 9 
photometric bands of the DIRBE instrument on COBE, and using the spec- 
trophotometry by the FIRAS instrument on COBE. 

The Japanese IRTS satellite obtained the 3-13 /im spectrum of the galactic 
plane {I « 45°, b « 0). If we assume that the 3-13 /im surface brightness is 
proportional to the 100 /im surface brightness, we can estimate the 3-13 /xm 
emission per H atom. The results are shown in Figure |S| The 3—13 /im 
emission appears to be dominated by the same infrared emission features 
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Fig. 8. Observed infrared emission per H nucleon from dust heated by the aver- 
age starlight background (see text). Crosses: IRAS (Boulanger & Perault 1988); 
Squares: COBE-FIRAS (Wright et al. 1991); Diamonds: COBE-DIRBE (Arendt 
et al. 1998); Heavy Curve: IRTS (Onaka et al. 1996, Tanaka et al. 1996). The 
interpolated dotted line is used to estimate the total power. 



observed in reflection nebulae and planetary nebulae, and widely attributed 
to polycyclic aromatic hydrocarbons. 

If we interpolate to obtain the dotted line in Figure|S| we can numerically 
integrate to obtain the total infrared emission power per H nucleon = 5.0 x 
10-2'*ergs-iH-i. About 2/3 of this power is radiated at wavelengths A > 
50 ^m, which is where interstellar dust at T ;S 20 K was expected to radiate. 
However, about 1 /3 of the power is radiated at A < 50 /im, where the emission 
from r < 20K dust should be negligible. Furthermore, about 1/6 of the 
total power is radiated at A < 10 /xm, primarily in 6.6, 7.7, and 8.6 /im PAH 
features. 



1.11 Interstellar Depletions: Atoms Missing from the Gas 

Narrow optical and ultraviolet absorption lines seen in absorption in stellar 
spectra can be used to determine interstellar gas phase abundances of many 
elements. While some species, such as N and S, have gas phase abundances 
(relative to H) which are approxiamtely solar, certain others, such as Mg, Al, 
Si, Ti, Ca, Fe, Ni, Cr show abundances which are far below solar. Since we 
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Fig. 9. Depletions in the diffuse molecular cloud toward (^Oph. 



presume that the interstellar abundances are approximately solar (Sofia & 
Meyer 2001), the atoms missing from the gas must be locked up in dust grains. 
In Fig. 121 the abundances have been plotted as a function of "condensation 
temperature" - the temperature at which cooling, solar abundance gas at 
LTE would begin to condense a particular element into a solid or liquid phase. 
Elements with Tcond ^ 1200 K are able to form "refractory" solids with large 
binding energies, and it is therefore not surprising that the strongly depleted 
elements tend to be those with high values of Tcond- The observed elemental 
depletions provide a clue to the composition of interstellar dust. Restricting 
ourselves to abundant elements, we see that the "electron donor" species in 
interstellar dust material must be predominantly C, Mg, Si, and Fe. The 
grains could contain a significant amount of O, and perhaps H. 

1.12 A Provisional Grain Model 

Below we will be discussing the astrophysics of interstellar dust. It is helpful 
to have a specific grain model in mind so that one can carry out calculations. 
Accordingly, here we describe briefly a grain model which appears to be 
consistent with observations. The model consists of a mixture of carbonaceous 
particles and silicate particles. The carbonaceous particles have the physical 
and optical properties of polycyclic aromatic hydrocarbon molecules when 
they are small a < 50A, or Nq ^ 6 x lO"* carbon atoms. When they are 




Fig. 10. WDOl size distribution for dust consistent with an Ry =3.1 extinction 
curve for the local Milky Way. The quantity plotted - a^dn/da - is proportional to 
the dust volume per logarithmic size interval. 



larger, the carbonaceous grains are assumed to have the optical properties 
of randomly-oriented graphite spheres. The silicate grains are assumed to be 
amorphous silicates. 

The emission in the 3—12 /j,m region from the diffuse ISM is attributed 
to dust grains small enough that a single starlight photon can raise the vi- 
brational temperature high enough for thermal emission in the observed vi- 
brational modes. In order for this emission to amount for ~ 21% of the total 
radiated power, very large numbers of very small dust grains are required. 
Li & Draine (2001) estimate that ~15% of the total carbon abundance must 
be in carbonaceous particle containing < 10^ C atoms. Since these particles 
must account for ~20% of the total absorption of starlight, it is obvious that 
they must make a signficant contribution to the interstellar extinction curve. 

Weingartner & Draine (2001a, hereafter WDOl) have shown that there are 
dust size distributions which include the required large numbers of ultrasmall 
carbonaceous grains and which are compatible with the observed extinction 
curves for various regions in the Milky Way, LMC, and SMC galaxies. The 
size distributions are by no means unique - the observed extinction curve 
cannot be "inverted" - but they are nevertheless strongly constrained, so the 
WDOl size distributions are probably not too far from the truth. Figure ITUl 



14 B. T. Draine 




2 4 6 8 

X-i (Atm-i) 

Fig. 11. Dashed line: average diffuse cloud extinction curve. Solid lines: model 
extinction curve for two extreme values of the ultrasmall carbonaceous grain abun- 
dance: be = (no PAHs) and he = 60ppm C in PAHs. For each case the size 
distribution of carbonaceous and silicate grains was adjusted to obtain good overall 
agreeement with the observed extinction. Prom Weingartner & Draine (2001a). 

shows a grain size distribution appropriate for diffuse clouds in the local Milky 
Way. The total mass in grains with a ^ lOOA is constrained by the observed 
extinction in the ultraviolet, but since these grains are in the "Rayleigh limit" 
the extinction is virtually independent of the actual sizes and numbers of 
these particles, provided only that the constraint on the total volume of 
these grains is satisfied. 

The size distribution for a ^ IsA carbonaceous grains is adjusted to make 
the predicted infrared emission agree with observations. The bimodal nature 
of the carbonaceous size distribution for a ^ lOOA is probably an artifact of 
the fitting procedure. 



1.13 Far-Infrared and Submm Opacities 

If dust grains in a cloud have a temperature Td, the emergent intensity is 

4 = B^Td) [1 - e-^'] , (2) 




Fig. 12. Opacities k(A) — (dust absorption cross section)/(mass of gas + dust) for 
Milky Way dust with Rv = 3.1 and Rv = 5.5, LMC average dust, and dust in 
the SMC bar, computed for the carbonaceous/silicate dust model (Weingartner & 
Draine 2001a; Li & Draine 2001). 



TdiX) = k(A) pds , (3) 



where the dust opacity 



p is the total mass density (gas + dust), and we have assumed riue/nu = 0.1. 
Dust opacities k are given in Fig. these opacities are also available at 
[http : //www . astro . princeton . edu/ ^draine/ dust/dust . html 

At long wavelengths the optical depth <C 1. If we can estimate the 
dust temperature, the cloud mass per area can be obtained from the observed 
intensity by 

p ds K — — . (5) 

k(A) B^iTd) 

The derived mass density depends sensitively on the value of Td unless one 
is in the Rayleigh-Jeans limit, A ^ hc/kTd, in which case J pds (x T^^ . 
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2 Optics of Interstellar Dust Grains 

Our knowledge of interstellar dust is based in large part on the interaction 
of dust grains with electromagnetic radiation: absorption, scattering, and 
emission. The interaction of electromagnetic radiation with a "target" will 
depend on the composition of the target (which determines the "response" of 
the target material to applied electric and magnetic fields) and the geometry 
of the target - both the shape and the size. Lacking definite knowledge of both 
the composition and the grain geometry, we must make some assumptions to 
proceed. 

2.1 Cross Sections, Scattering Matrices, and Efficiency Factors 

Wc arc generally interested in calculating the following quantities for unpo- 
larizcd plane waves incident on the target; 

• C'abs = the total absorption cross section. 

• C'sca = the total scattering cross section. 

• dCsca/d^ — differential scattering cross section. This is related to the 
MuUer matrix element by dCsca/d^ = Sn/k"^, where k = 2tt/X. 

• polarization P for the radiation scattered in a particular direction. 

In some cases we wish to explicitly consider polarized light. For this case 
the Stokes vector for the scattered radiation is obtained by multiplying the 
Stokes vector of the incident radiation by the 4x4 Mueller scattering ma- 
trix Sij. The 4x4 scattering matrix (which, for a fixed grain, is a function of 
direction of incidence and direction of scattering) fully characterizes the scat- 
tering properties of a grain. Sec Bohren & Huffman (1983) or Mishchencko, 
Hovenier, & Travis (2000) for discussions of scattering concepts. 

It is convenient to normalize scattering and absorption cross sections by 
dividing them by a geometric cross section. In the case of spherical target 
of radius a, the standard convention is to simply divide by the geometric 
cross section ira^. In the case of nonspherical targets, there is more than one 
convention in use. Some authors divide by the geometric cross section of the 
target as seen from the direction of the illumination. Other authors divide 
by the average geometric cross section for random orientations. 

My preference is to normalize by dividing by the geometric cross section 
of an "equal volume sphere" . Thus for a target of mass M, and made of a 
material with density p, the volume is y = M/p, and we define "eflaciency 

factors Qsca) Oabs) Qext = Osca "I" Oabs by 



Q = 



c 



'eff 



where 




(6) 
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2.2 Grain Geometry? 

Since interstellar dust is able to polarize starlight, we know that at least some 
interstellar dust grains must be appreciably nonspherical. Nevertheless, for 
many purposes we will approximate dust grains as spheres for the simple 
reason that we can use "Mie theory" to calculate scattering and absorp- 
tion of light by spherically-symmetric targets with sizes comparable to the 
wavelength. Exact series solutions have been found for spheroids (Asano & 
Yamamoto 1975, Asano & Sato 1980; Voshchinnikov & Farafonov 1993) but 
numerical calculations are much more challenging than for spheres. Solutions 
for other shapes (other than the unphysical case of infinite cylinders) are not 
available. 



2.3 Dielectric Functions 

First, a few words about conventions: We will follow c.g.s. electromagnetism, 
so that e = /i = 1 for the vacuum. We will represent plane waves as propor- 
tional to g^kx-iujt.^ with this convention, absorption corresponds to Im(e) > 0. 

When an electric field is applied to a material, there are two distinct types 
of response possible: the response of the bound charge and the free charge. 
The bound charge (e.g., electrons and nucleus within an atom) undergoes 
a finite displacement when the electric field is applied, with the individual 
atoms or molecules acquiring an electric dipole moment in response to the 
applied electric field. The response of the bound charge is characterized by a 
complex dielectric function e(''°""'')(aj). 

The free charge responds in the form of an electric current density, J = 
aE, where ct{oj) is the electrical conductivity. 

When dealing with a monochromatic electric field it is convenient to define 
a dielectric function e describing the response of both bound and free charge: 

+ . (7) 

LU 

With this definition, the free charge current J is now absorbed into the dis- 
placement current {l/4TT)dD/dt, so that Maxwell's equation becomes 

where D — eE. 

[Oi l A Model Insulator 

It is instructive to consider a simple classical model for the dielectric 
response of an insulator. Suppose the insulator consists of "molecules" , with 
no permanent dipole moment, but which can be polarized in response to 
the local applied electric field plus the electric field from all of the other 
polarized molecules. The dipole moment of a molecule is qx, where x is the 
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displacement of the bound charge q. Suppose the displacement x behaves like 
a driven, damped harmonic oscillator, 

mx ^ qE — mujQX — mx/ro , (9) 

where is the resonant frequency, and tq is a damping time. For a periodic 
electric field E oc e~''^* we can readily solve for the molecular polarizability 
a = qx/E: 

"(t^) = -; r~/ — ■ If 2 \ ' (10) 

where ^ 

ao = (11) 

is the zero- frequency polarizability of the molecule. We now consider a medium 
consisting of a number density n of such molecules. Each now responds to 
the externally applied electric field plus the electric field due to all of the 
other polarized molecules. In the low frequency limit, the effective dielectric 
function e of the medium is related to the molecular polarizability through 
the famous Clausius-Mossotti relation (see, e.g., Jackson 1975) 

Anna , 
1 — (47r/3)nQ! 

For molecules on a cubic lattice, the Clausius-Mossotti relation (|12|l is exact 
in the limit |m|A:(i <C 1, where d = n~^/'^ is the intermolecular spacing and 
k = oj/c. For finite kd, corrections to 0[{kdf'] have been derived by Draine 
& Goodman (1993). Substituting (fTHI into itT^ we obtain 



1 - (47r/3)nao - {uj/looY - «^^/(^o'^o) 



In the low frequency limit — s- 0, we have 



Re(e(''°""'^))^4^°""'^)^l+ , (14) 

1 — (47r/3)Q;o 



from which we note that 

47rao =: 3 



{hound) -| 

fo 

{hound) . rj 



(15) 



At low frequencies, the imaginary part of e(''°""'') varies linearly with fre- 
quency: 

jj^(g(6o«„d)) « Aoj , (16) 

Q {bound) -| 
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An insulator has no mobile charges (tr ~ 0), and for our model at low fre- 
quencies we have 

Re(e) ^e^"^ =consi , (18) 
Im(e)^Acj . (19) 

I2T3I 2 A Model Conductor 

For a material with a density rig of free electrons, a simple classical model 
for the electron dynamics would be 

rUeX ~ eE — m^x/Te , (20) 

where is the electron collision time. With this equation of motion, the 
conductivity (7{lu) = rieex/ E is just 

a(c.) = -^^ , (21) 

1 — tUJTe 

where 

do = Te/rrie = — — , (22) 

where ctq is the d.c. conductivity, and ujp = (A-Knee^ /me)^^^ is the plasma 
frequency. For this case, the low-frequency behavior of e is 

Re (e) — + ^(bound) _ ^^^^^^ _ ^gj^gf ^ ^23) 

, s . 47r(Tn 

Im (e) Alu + . (24) 

Comparing p9|) and (|24l) , we see that the low- frequency behavior of Ime 
is qualitatively different for an insulator or for a conductor. Below we will see 
what this implies for the frequency-dependence of the absorption by small 
particles. 



2.4 Calculational Techniques 

In many astrophysical applications we are concerned with particles which are 
neither very small nor very large compared to the wavelength of the incident 
radiation. There are several different methods which can be used to calculate 
scattering and absorption cross-sections by targets which are neither very 
large nor very small compared to the wavelength: 

• Mie theory solution for spherical targets (Bohrcn & Huffman 1983). Mie 
theory codes are highly-developed and readily available. They break down 
due to numerical roundoff error when the target becomes too large com- 
pared to the wavelength. A modified version of the Bohren & Huffman 
code is available at 

ihttp : / / www ■ astro . princeton . edu/ ^draine/ scattering . html^ 
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• Series solution for homogeneous spheroid or layered spheroid (boundaries 
must be confocal spheroids) (Asano & Yamamoto 1975; Asano & Sato 
1980; Voshchinnikov & Farafonov 1993). Codes for spheroids have been 
developed but have not been widely used. Codes based on the Asano 
& Yamamoto treatment had a reputation for being somewhat numeri- 
cally delicate (Rogers & Martin 1979). Voshchinnikov & Farafanov (1993) 
have developed a treatment based upon separation of variables which is 
reported to be more robust. They have recently generalized this to core- 
mantle spheroidal targets (Farafanov, Voshchinnikov, & Somsikov 1996). 

• Extended Boundary Condition Method (EBCM), often referred to as the 
"T-matrix method" (Mishchenko, Travis & Macke 2000). The EBCM 
method (Mischchenko, Travis & Macke 2000) is used to construct the "T 
matrix" which gives the coupling between vector spherical harmonic com- 
ponents of the incoming wave and vector spherical harmonic components 
of the scattered wave. Once the T-matrix has been obtained, orienta- 
tional averages can be efficiently calculated. The EBCM method appears 
to be well-suited to targets with rotational symmetry which are not ex- 
tremely elongated. EBCM T-matrix codes have been made available by 
M. Mishchenko at jhttp : //ww w . giss . nasa . gov/ ^crmim/tjn atrix . html] 

• Discrete Dipole Approximation (Draine & Flatau 1994; Draine 2000). The 
Discrete Dipole Approximation can be quite readily applied to complex 
geometries. It has been used to calculate the optical properties of graphite 
particles, with a highly anisotropic dielectric tensor (Draine 1988; Draine 
& Malhotra 1993). In recent years it has been accelerated by use of EFT 
techniques (Goodman, Draine, & Flatau 1991), and its accuracy has 
been improved by refinements in the assignment of dipole polarizabili- 
ties (Draine & Goodman 1993). The Fortran code DDSCAT is available 
from | http : //www . astro . princeton . edu/ ^draine/DDSCAT . html , and 
a detailed User Guide is available (Draine & Flatau 2000). 

2.5 Scattering by Homogeneous Isotropic Spheres 

The solution to Maxwell's equations for a plane wave incident on a homo- 
geneous and isotropic sphere was independently obtained by Mie (1908) and 
Debye (1909). The solution - now generally referred to as "Mie theory" - 
is given in terms of a series expansion in powers of the "size parameter" 
X = 27ra/A. The expansion converges, but the number of terms which must 
be retained is of 0{x). The details of the Mie theory solution are nicely de- 
scribed in the excellent monograph by Bohren & Huffman (1983). There are 
a number of computer codes to evaluate the absorption cross section and 
differential scattering cross sections given by Mie scattering theory. When x 
is large, care must be taken to avoid errors due to finite precision arithmetic 
in numerical evaluation of the series solution. 

For an incident monochromatic plane wave, the full Mie theory solution 
depends on the dielectric function e(u}) - describing the electric polarization 



Astrophysics of Dust in Cold Clouds 



21 



of the material in response to an applied electric field oscillating at angular 

frequency uj and the magnetic permeability function /i(cL)) characteriz- 
ing the magnetization of the material in response to an applied oscillating 
magnetic fields. Because the magnetic response is generally negligible (i.e., 
l/i — 1| <C |e — 1|) at frequencies above 100 GHz) it is customary to neglect 
the magnetization of the grain material except insofar as it is due to eddy 
currents. 



y 




tetrahedron 
N=61432 

Fig. 13. Representation of a tetrahedron with a dipole array. The dipoles are actu- 
ally pointlike - here they arc represented by small spheres for purposes of visualiza- 
tion. We consider incident radiation propagating along the direction of the vector 
kin illustrated here. Taken from Draine (2000). 



2.6 Discrete Dipole Approximation 

The discrete dipole approximation consists in replacing the (continuum) tar- 
get of specified dielectric function e with an array of polarizable points (re- 
ferred to as "dipoles"). The array geometry is chosen to mimic the target 
geometry. The polarizabilities a of the points are chosen so that an infinite 
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Fig. 14. Sii = k^dCeca/dfl for scattering of unpolarized light bya tetrahedron, 
for scattering directions in the x-y plane (see Fig. 113^ , as a function of scattering 
angle O. The tetrahedron has refractive index m = 1.33 + O.Oli, has edges of length 
s — 7.849A, where A is the wavelength of the incident light in vacuo, and is oriented 
with an angle 9 — 30° between kin and ai. The peak at © = 240° corresponds 
to the direction of specular reflection for geometric optics. Results are shown for 
a tetrahedron represented by A*' = 61432 and N = 105832 dipoles. The scattered 
intensities are in excellent agreement in all directions where the scattering is at all 
strong. The dashed line shows k^{Zii) = Csca/47r. Taken from Draine (2000). 

lattice of such polarizable points would have the same dispersion relation as 
the material of dielectric function e (Draine & Goodman 1993). With FFT 
techniques employed to speed the calculation, it is now feasible to calculate 
scattering and absorption by targets represented by more than 100,000 po- 
larizable dipoles on a workstation with 256 MB of RAM. There are great 
efficiencies if the dipoles are situated on a cubic lattice, so DDSCAT requires 
that this be the case. 

The DDA can be applied to inhomogeneous targets and targets with com- 
plex geometries. As an example of the type of problem which can be solved, 
in Fig. ^|we show a discrete-dipole array of 61432 dipoles intended to ap- 
proximate a tetrahedral target. In Fig.^jwe show the calculated scattered 
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intensity in the x-y plane as a function of scattering angle Q. Notice the pro- 
nounced scattering peak at = 240° - this is the direction where one would 
have a specular reflection peak in the geometric optics limit. 

For a fixed target size D, complex refractive index m, and wavelength A, 
the DDA converges to the exact answer in the limit where the interdipole sep- 
aration d ^ 0: the target structure is well-resolved and the dipole separation 
is small compared to the wavelength in the target material. The criterion for 
separation small compared to A can be written \m\kd ^ 1. If |m| is not too 
large, we find that the overall scattering and absorption cross section is cal- 
culated fairly accurately provided \m\kd ^ 1; the detailed scattering pattern 
will be calculated accurately if Imj/cd ^ 0.5. For example, the two calculations 
in Fig. El with Imj/crf = 0.82 and 0.68, are already in fairly good agreement. 



2.7 Infrared and Far-Infrared 

When the wavelength is long compared to the target size, we may use the 
"dipole approximation" (Draine & Lee 1984) and characterize the response 
of the target to the incident electromagnetic field purely in terms of the 
(oscillating) electric and magnetic dipole moments of the target. We may 
write 

Cabs = — [Im(«|) + Ini(aS)] , (25) 
c 

Csca=f (^^)'[KP + |aSn , (26) 

where a|; is the (complex) electric polarizability along the direction of the 
incident electric field vector, and is the (complex) magnetic polarizability 
along the direction of the incident magnetic field vector. 

Consider a homogeneous ellipsoidal particle composed of an isotropic ma- 
terial, with semiaxes a,b,c. If the applied electric field is along one of its 
principal axes, the electric polarizability is 



V £-1 

4^(e-l)L£; + l 



(27) 



where Lb is the "shape factor" along the direction of the applied E field, and 
V = Anabc/S is the volume of the ellipsoid. The shape factors along the three 
principal axes satisfy the sum rule Li + L2 + L3 — 1 and may be obtained 
by numerical quadrature (Bohren & Huffman 1983). 

Consider now a spheroid with semiaxes a, b, b (prolate if a/5 > 1, oblate if 
a/b < 1). For this case we have analytic expressions for the shape factors La 
and Lb = (1 — La)/2. For a prolate spheroid (a < b) we have (van de Hulst 
1957) 

La = 



^lnfi±^)-l 
2e V 1 - 



(28) 
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and for an oblate spheroid (a < b) we have 



L 

where 



1 

1 arctan(e) 

e 



(29) 



e'^\l-ib/ar\ . (30) 

Exact solutions for a% in the limit A 3> Ocfr are also available for layered 
grains provided the interfaces are confocal spheroids (Gilra 1972; Draine & 
Lee 1984). 

For spheres we have La = = 1/3, and 

^"•^^^ c (ei+2)2 + e2 ' (31) 

where ei = Re(e) and £2 = Ini(e). 

Eq. (|31|l contains an important result: the absorption cross section of a 
grain with size ^ A is simply proportional to the grain volume! The far- 
infrared opacity therefore depends only on the total volume of grain material 
present, but not on the sizes of the particles, provided only that they are 
small compared to the wavelength A. 

How do we expect the absorption cross section to depend on frequency uj 
at low frequencies? For an insulator, from 119() we see that 

W A 

while for a conductor, we see from (|24|l that 

Cabs- — f-^) . (33) 

Thus we see that for both insulating grains and conducting grains, a simple 
physical model leads to Cabs « uo"^ at low frequencies! This is the basis for 
the expectation that dust opacities should vary as A^^ in the far-infrared, as 
in the dust models of Draine & Lee (1984). 

Now the above discussion has been based on simple classical models for 
the response of charge to applied electric fields. Real materials obey quantum 
mechanics, so results may differ. Furthermore, the above discussion assumed 
that w <C cjo in discussion of the response of the bound charge, but an 
amorphous material could have very low frequency vibrational modes, so 
that the assumption that w <C wo might not be valid in the far-infrared. 
This might explain why some laboratory studies of amorphous materials find 
different behavior: for example, Agladze et al. (1996) find that the absorption 
coefficient for amorphous MgO-2Si02 varyies as w^'^ between 800 /im and 4 
mm. 

Nature must be our guide, but it is worth keeping in mind that opacities 
varying as uJ^ emerge naturally from simple models. 
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2.8 Kramers-Kronig Relations 

PurccU (1969) pointed out that the Kramers-Kronig relations can provide 
useful constraints in dust modelling. The Kramers-Kronig relations are gen- 
eral relations which apply to "linear response functions" , such as a dielectric 
function, which specifies the response (e.g., the electric polarization) to an 
applied stress (e.g., an applied electric field). The only assumptions are that 
(1) the response is linear, and (2) the system is causal - the response can 
depend on the stress applied in the past, but cannot depend upon the future. 
With these very simple assumptions, it is possible to derive the Kramers- 
Kronig relations (see Landau & Lifshitz 1960 for a nice derivation). In the 
case of the dielectric function, the Kramers-Kronig relations are 

eiiiOo) = l + -P r duj^^^^ , (34) 
e^iiuo) ^ -iOoP r diu^^ , (35) 

where P indicates that the principal value is to be taken. 

Thus the real and imaginary parts of e{uj) are by no means independent 
- if either one is known at all frequencies, the other is fully determined. In 
order to have a physically acceptable dielectric function, one approach (see, 
e.g., Draine & Lee 1984) is to specify £2 at all frequencies, and then construct 
ei(tj) using 



2.9 Kramers-Kronig Relations for the ISM 

Purcell showed that one could apply the Kramers-Kronig relations directly 
to the interstellar medium. Plane waves propagate through the interstellar 
medium but are attenuated by scattering and absorption. One can describe 
this attenuation by an imaginary component of the dielectric constant e of 
the interstellar medium, where |e — 1| <C 1. 

Recall that for a electromagnetic plane wave of frequency the wave 
vector is k = m{uj)u;/c, where m{uj) — ^/e is the complex refractive index. 
The electric field intensity decays as 

E cx exp [-Im(fc)x] = exp [-(l/2)Im(g)w/c] . (36) 

The energy in the wave is proportional to \E\'^ and therefore decays twice as 
rapidly as E. Thus 

ngrCcKtW = —£2 ■ (37) 

We can now apply the Kramers-Kronig relation (|34|l to obtain the real 
part of the dielectric function e of the ISM at zero frequency: 

ei(0)-l = -/ — -%.Cext(w) . (38) 




Fig. 15. Function F{a/b,eg) as a function of the static dielectric constant eo, for 
selected values of axial ratio a/b. After Purcell (1969). 



The static dielectric function of the ISM is directly related to the electric 
polarizability agr of the grains in it: 



il{0) — 1 ~ ATTngragr{0) 



Thus 



4:Trngragr{0) = —cug 



1 



(39) 

(40) 
(41) 



We have seen above how the static polarizability of a spheroidal grain de- 
pends on its volume, shape, and dielectric function. Averaged over random 
orientation, we have (for a dielectric grain) 



47ragr(0) = 31/F(a/6,eo) , 



F{a/b, eo) 



(eo-l)3L„ + 3 (eo- 1)3^6 + 3 



(42) 
(43) 



F is just the orientationally-averaged polarizability relative to the polariz- 
ability of an equal- volume conducting sphere. Using H42(l . we obtain the grain 
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volume per H atom in terms of an integral over the extinction per H atom: 
^-t4t; rdA^Ce.t(A) (44) 

Now consider a conducting material. From (|23I24|I we see that e ix oo, 
so that 

so that a conducting sphere {La — Lb ~ 1/3) has F{a/b = 1, ep = oo) — 1- 

We do not know the grain shape, and we do not know the static dielectric 
constant eg. However, from Figure [T31 we see that F < 1.5 unless the grain is 
extremely elongated {a/b > 20 or b/a > 20) and the static dielectric function 
is very large (eo > 10) - as for a metal. For a reasonable grain shape (e.g., 
0.5 < a/b < 2) and dielectric function cq ~ 3 we have F « 0.4. 

Purcell's analysis is a delight, and has two important consequences: 

12.91 1 Grain volume per H 

It is of course not possible to measure the extinction per H atom at 
wavelengths from to oo. However, since Cext > 0, measurements over a 
finite wavelength range can be used to obtain a lower bound on FVugr/n-n. 
The extinction per H nucleon is fairly well-known from 0.1 /^m to 30 /im; a 
numerical evaluation gives 



/ Z^^ii^dA^ 1.1 X 10-25 cmVH . (46) 

J 0.1 Mm 



Approximately 50% of the integral is contributed by wavelengths 0.1 < A < 
1 /xm. If we estimate F k, 0.4 this gives 

^ > 9.3 X 10-2^ cmVH , (47) 
riH 

or, for an assumed grain mass density p = 2.5 gcm^^ (intermediate between 
graphite and silicate) we have a lower bound on the ratio of grain mass to H 
mass: 

^>.014 . (48) 

I2.9L 2 Asymptotic Behavior of Ccxt at Long Wavelengths 

Eq. (|44|l tells us that Ccxt(A)(iA must be convergent, and therefore that 
Ccxt (A) must decline more rapidly than 1 /A: grain models in which Cext oc 1/A 
for A — > oo are unphysical. It is of course possible to have d In Ccxt /^^ In A ~ 
— 1 over a limited range of wavelengths, but this cannot be the asymptotic 
behavior. 
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Fig. 16. Microwave sky brightness correlated with 100 nm emission. Shown is the 
estimated microwave emissivity per unit H nucleon. 



2.10 Microwave 

Dust grains would be expected to radiate thermally at microwave frequencies, 
but this emission was expected to be quite weak, based on an extrapolation 
from the thermal emission peak at ~ 100 fim. It was therefore a surprise 
when sensitive maps of the microwave sky brightness (for the purpose of 
studying angular structure in the cosmic background radiation) showed rel- 
atively strong emission correlated with the Galactic 100 fim emission, and 
therefore with interstellar dust (Kogut et al. 1996a, b). Figure [TBI shows the 
inferred microwave emissivity of per H nucleon. The dashed line labelled 
"/3 = 1.7, Td = 19.5" is the contribution expected from thermal emission 
from dust grains if the dust opacity oc A""'^''' for A ^ 100 /Ltm. We see that at 
30 GHz (A = 1 cm the observed microwave emission is two orders of magni- 
tude above the value expected from classical dust grains. 

The power radiated in the microwave is a tiny fraction of the total, so this 
microwave emission is not an important source of cooling. Nevertheless, we 
would like to know what process is responsible. Possible mechanisms include 
synchrotron emission from relativistic electrons, or free-free emission from a 
thermal plasma, but these seem unable to account for the observed 15-50 
GHz emission (see Draine & Lazarian 1999b). We are left with interstellar 
dust. 



Astrophysics of Dust in Cold Clouds 



29 



The possibility of radio emission from rotating dust grains dates back 
to Erickson (1957). As discussed above, the strong infrared emission in the 
3 — 12 /xm range appears to require a very large population of very small dust 
grains. Ferrara & Dettmar (1994) pointed out that these dust grains would 
produce detectable radio emission if they were undergoing Brownian rotation 
at the gas temperature. 

Draine & Lazarian ( 1998a, b) analyzed the grain rotational dynamics, and 
showed that the population of ultrasmall grains required to understand the 
infrared emission would be expected to produce microwave emission, with 
a predicted spectrum (and intensity) as shown in Fig. El We will discuss 
the grain dynamics in ^ but this rotational emission still seems likely to 
contribute a substantial fraction of the observed dust-correlated microwave 
emission. In principle, these rapidly spinning dust grains could be aligned, 
leading to polarized emission, but the polarization fraction has been esti- 
mated to be small (Lazarian & Draine 2000). 

However, if dust grains contain magnetic materials - which does not seem 
implausible, considering the large amount of Fe in interstellar grains - there 
could also be thermal magnetic dipole emission from dust. Consider a ferro- 
magnetic domain in a dust grain. This domain is spontaneously-magnetized: 
it is energetically favorable to have alignment of electron orbital angular mo- 
mentum and electron spins. The lowest energy state has some specific mag- 
netization, but there are nearby energy states where the domain has uniform 
magnetization with the same net magnetic moment but in a slightly differ- 
ent direction. In thermal equilibrium, it is possible for these higher energy 
states to be excited, so that the magnetization direction will fluctuate. This 
results in a time- varying magnetic dipole moment, and hence magnetic dipole 
radiation. 

The optics of magnetic dust grains with size small compared to the wave- 
length has been developed by Draine & Lazarian (1999a; hereafter DL99a), 
and the magnetic susceptibility /x(w) was estimated for several magnetic ma- 
terials, including metallic Fe and magnetite Fe304. The spectrum and in- 
tensity of this radiation was estimated by DL99a, and it was found that 
if a substantial fraction of the Fe in interstellar grains is incorporated into 
magnetic materials (e.g., magnetite), the resulting thermal emission could 
account for an appreciable fraction of the observed microwave emission! 

How will we be able to distinguish between spinning dust grains and 
magnetic dust grains as sources of microwave emission? DL99a point out 
that small grains are believed to be underabundant in dense clouds (they 
presumably coagulate to form bigger grains), so that the rotational emissivity 
per H nucleon would be reduced. The magnetic dipole emission from magnetic 
grain materials does not, however, depend on the grain size. Hence if we see 
strong microwave emission from dust in dense clouds, this would suggest 
that magnetic dust grains may be responsible. Future pointed observations 
of dense clouds may answer this question. 
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2.11 X-rays 

The scattering of X-rays by dust grains was first discussed by Overbeck 
(1965), Slysh (1969), and Hayakawa (1970). Because the refractive index of 
grain materials is close to unity at X-ray energies, the scattering is through 
small angles. As a result, images of X-ray point sources show a "halo" of 
dust-scattcrcd X-rays. The observable halo can extend for tens of arcminutes 
from the source. 

For sufficiently small grains, this scattering can be calculated using the 
Raylcigh-Gans approximation (sec Bohren & Huffman 1983), which is valid 
if \m — l\ka <C 1, where a is the grain radius. Because X-rays have large 
k, this condition may not be satisfied for the larger interstellar grains, in 
which case one should resort to full Mic theory. This can be numerically 
challenging because the number of terms which must be retained is of or- 
der ka ~ 507(a/0.1 /um)(£'/keV), so high numerical accuracy is required. 
Smith & Dwck (1998) have compared scattering halos estimated using the 
Rayleigh-Gans approximation with Mie theory, and show that the Rayleigh- 
Gans approximation fails for energies below 1 keV. Scattering by grains with 
|m — l|/ca <C 1 and ka ^ 1 can be treated using "anomalous diffraction 
theory" (van de Hulst 1957). 

Nova Cygni 1992, a bright X-ray nova, was observed by ROSAT. The 
observed X-ray halo has been compared to predictions for different grain 
models by Mathis et al. (1995), Smith & Dwek (1998), and Witt, Smith & 
Dwek (2001), with differing conclusions. The most recent study (Draine & 
Tan 2002) concludes that our standard dust model is in reasonable agreement 
with the observed X-ray halo. Future observations of X-ray scattering halos 
will be valuable to test and constrain grain models. 



3 IR and Far-IR Emission from Interstellar Dust 
3.1 Heating of Interstellar Dust 

While dust grains can be heated by collisions with gas atoms and molecules, 

starlight usually is the dominant heating process. The rate of energy depo- 
sition, the rate of photon absorptions, and the mean energy per absorbed 
photon are 

/•OO 

{dE/dt)ehs = dv c Cabs(z^) , (49) 
Jo 

Nabs = J dv-^cCabs{v) , (50) 

{dE/dt)abs ,^y. 

(Mabs = , (51) 

jVabs 

where Uvdv is the photon energy density in {u, u+du), and Cabs{f) is the pho- 
toabsorption cross section. We assume that all of the energy of an absorbed 
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photon is converted to heat; this is not exactly correct, since energetic pho- 
tons may eject a photoelectron, or excite fluorescence, but these processes 
take away only a small fraction of the total absorbed power. 

Now suppose that the dust grain has energy E distributed among its 
vibrational degrees of freedom. Except when the energy E is very small, the 
number of vibrational modes which can be excited is very large. The density 
of states has been discussed by Draine & Li (2001): even for a small molecule 
like coronene C24H12, there are ~ 10^° different states with total energy 
E < leV; for C4000H1000 this number becomes ~ 10^^^. As a result, the 
statistical notion of "temperature" can be used, even for quite small grains: 
a grain with vibrational energy content E is assumed to be characterized by 
a temperature T{E) such that if the grain were in contact with a heat bath 
at temperature T, the expectation value for its energy would be E. 

At temperature T, the average power radiated by a grain is 



where Cabs('^) is the angle-averaged absorption cross section for photons of 
frequency and 



is the Planck function. 

It is now natural to determine the "steady-state" temperature Tgs at which 
the power radiated equals the power absorbed: 



Once we solve for Tgg, we calculate the thermal energy content of the grain 
at this temperature, E(Tss)- There are two regimes: 

• If E{Tss) S> (/ii^)abs: then individual photon absorption events do not 
substantially change the energy content of the grain, and we may assume 
that the grain temperature T = Tgs- This is the "steady heating" regime, 
where we can approximate the discrete heating events as a continuous 
process. 

• If E(Tss) <C {hi')ahs, then individual photon absorption events heat the 
grain up to peak temperatures Tmax ^ Tss, and the grain will usually cool 
to a temperature Tmin < Tss before the next photon absorption event. 
This is the "stochastic heating" regime, where the discrete nature of the 
heating has important consequences. As we will see below, stochastic 
heating is often important. 




(52) 




(54) 



Figure [T7I shows the temperature history of 4 grains of different sizes in the 
diffuse interstellar medium over a ~ 10^ s interval - about one day. The 
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Fig. 17. A day in the life of an interstellar grain: grain temperature vs. time for 
4 grain sizes, for grains heated by the average interstellar radiation field. Grains 
with a 200A have a nearly constant temperature, but a ^ lOOA grains show con- 
spicuous increases in temperature following each photon absorption, with gradual 
cooling between photon absorption events. 



temperature of the a = 200A grain fluctuates in a small range around ~ 20 K. 
However, as the grain size is reduced, the fluctuations become increasingly 
extreme, with the peaks being higher and the low points being lower. The 
smallest grain shown (a = 25A) absorbs about 1 starlight photon per day; this 
grain spends most of its time quite cool, but immediately following absorption 
of a photon the grain temperature reaches ~ 50 K. The grain cools by infrared 
emission; it is obvious that most of the infrared emission must take place 
during the brief interval while the grain is "hot" . If we wish to calculate 
the time-averaged emission spectrum, we cannot use the "average" grain 
temperature - we need to integrate over a distribution of grain temperatures. 
This phenomenon is even more pronounced for smaller grains. A grain with 
a radius of SA has a volume only 1/125 that of the 25A shown in Figure [T7I 
the mean time between photon absorptions is of order ~ 10^ s, but a 10 eV 
starlight photon can raise the grain temperature to ~ 10"^ K. 
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Let P{E) be the probability of the grain having vibrational energy E' > 
E. In the steady heating regime, we approximate P{E) as a step function, 
and \dP/dE\ as a delta function 5{E — Egg)- 




Fig. 18. Energy distribution functions for charged carbonaceous grains with radii 
a = 5, 10, 25, 50, 75, 100, 150, 200, 250, 300A in the interstellar radiation field. The 
discontinuity in the 5, 10, and 25A curves is an artifact due to a change in the 
method of estimating the cooling when the energy is equal to the 20th vibrational 
mode. For 5, 10, and 25A, a dot indicates the first excited state, and the probability 
Pq of being in the vibrational ground state is given. Taken from Li & Draine (2001). 



In the stochastic heating regime, we must solve for P{E). Draine & Li 
(2001a) show examples of energy distribution functions. If the energy states 
arc grouped into bins j = 0,1,... (where j = is the ground state), then 
we can calculate the state-to-state transition rates Tjj for transitions i ^ j 
due to both photon absorptions and photon emissions. Once Ti^j is known, 
we define the diagonal elements Tu = ^^j^^Tji. The steady state solution 
Pj for the probability of being in state j then satisfies the N coupled linear 
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equations 

M 

= ^ T.jPj for i = 0, M . (55) 

3=0 

Using the normalization condition X^^lo -^j ^ ^ ' obtain a set of M linear 
equation for the first M elements of Pj: 

M-l 

^{T^J-T^M)P3^~TaI for i = 0, . . . , M - 1 , (56) 
i=o 

which we solve using standard techniques. In practice, we take M w 500, in 
which case iterative methods are required to efhciently solve for Pj . Once the 
Pj are determined for j — 0, ...,M — 1 we could obtain Pm by subtraction 
Pm = 1 — Sjf=o^ t^^^ inaccurate; it is better to instead use 

j=o 

Fig. ^] shows dP/d In E for carbonaceous grains heated by the average 
intestellar radiation field (ISRF) due to starlight in an H I region, as esti- 
mated by Mathis, Mezger, & Panagia (1983, hereafter MMP). For a < 25A 
grains dP/dlnE becomes very small for E/hc > (l/91lA) = 1.1 x lO'^cm^^: 
this reflects the fact that grains cool essentially completely between photon 
absorption events, so that the energy content virtually never rises above the 
maximum energy (13.6 eV) of the illuminating photons. As the grain size 
increases, the time between photoabsorptions goes down, and the cooling 
time at fixed energy goes up; when the grain size exceeds ~ 30A there is a 
significant probability of a photoabsorption taking place before the energy 
of the previous photoabsorption has been radiated away. However, stochastic 
heating effects are noticeable even for grains as large as 300A - the energy dis- 
tribution function has narrowed considerably, but is still appreciably broad. 



3.2 IR and Far-IR Emission Spectrum 

With the energy distribution function calculated as discussed above, we can 
now calculate the time-averaged emission spectrum for a carbonaceous grain: 

f dP 

F^^in J dE—C,U>^)B^{TiE)) . (58) 

In Fig. 1191 we show the emission spectrum of PAII+ molecules of various sizes 
heated by the ISRF, and in Fig. [201 we compare the emission spectrum calcu- 
lated sing the energy distribution functions dP/d In E with emission spectra 
calculated for dust with steady temperature T^s- We see that stochastic heat- 
ing is important even for grains as large as ~ lOOA in the ISRF. 
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X (/um) 

Fig. 19. IR emissivity per C atom for PAH"'" molecules of various sizes illuminated 
by the average ISRF. Solid line: "exact-statistical" calculation, which uses tran- 
sition probabilities which do not involve the thermal approximation. Dotted line: 
"thermal-discrete" calculation where the spontaneous emission rates are calculated 
using a thermal approximation, as discussed in the text. The two methods are 
in excellent agreement, indicating that the thermal approximation can be used to 
calculate the transition rates. Taken from Draine & Li (2001). 

In Fig. 1211 we show the emission spectrum calculated for a mixture of 
carbonaceous and silicate grains with size distributions consistent with the 
average (Ry = 3.1) interstellar extinction curve, illuminated by the local 
ISRF. This should be the emission spectrum of the "cirrus" clouds. Also 
shown in Fig. 1211 are observations of the emissivity per H atom of the dust 
at high galactic latitudes. 

Diamonds show DIRBE photometry; the small triangles show our calcu- 
lated spectrum convolved with the DIRBE filters, for comparison with the 
diamonds. The squares show the FIRAS determination of the emissivity per 
H nucleon. The far-infrared emission is in excellent agreement with the pre- 
dictions of our dust model. The model is in very good agreement with DIRBE 
photometry at 3.5 /im and 25 fiui; the model appears too low by about a fac- 
tor 2 at 5 and 12 /im. The observations are difficult; it is hoped that SIRTF 
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Fig. 20. IR emission per grain for carbonaceous grains of various sizes illuminated 
by the average ISRF. Solid line: full stochastic heating calculation. Broken line: 
assuming steady-state temperature Tss- For a ^ 200 A in the ISRF we see that 
stochastic heating has little effect on the emission spectrum, but for a ^ 150A it 
significantly modifies the emission spectrum. Taken from Li & Draine (2001). 

will be able to measure the spectrum of selected cirrus clouds for comparison 
with our model. 

In Fig. 1221 we compare our model with observations taken on the galac- 
tic plane, where the FIR surface brightness is much higher. In addition to 
DIRBE photometry, spectroscopic observations made with the IRTS satellite 
are available. 

These observations are looking in the galactic plane at the "tangent" point 
of gas in the 6 kpc ring, at a distance of 6 kpc from us. At this distance, the 
40' DIRBE beam subtends 70 pc, so the observations sum the emission from 
dust in a region ^70x70 pc on a side, and many kpc long. We compare the 
observations to a model which is obviously oversimplified: a uniform slab 
of dust heated by a uniform radiation field with the spectrum of the local 
ISRF but whose intensity we scale so that the dust emission best reproduces 
the observed spectrum: the best fit is obtained with the starlight intensity 
equal to twice the ISRF. The DIRBE beam obviously includes emission from 
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Fig. 21. IR emission per H for dust mixture illuminated by the ISRF. Diamonds: 
DIRBE (Arendt et al 1998). Squares: FIRAS (Finkbeiner et al. 1999). From Li & 
Draine (2001). 

stars, which we believe dominate at A < 4 /xm, so our model includes light 
from stars assumed to be mixed uniformly with the dust. This simple model 
is in good agreement with the observations. In particular, we successfully 
reproduce the very strong emission in the 5-12 ^m region. The discrepancy at 
100 fiui is somewhat surprising, but may in part be due to oversimplification 
of our model, and perhaps in part to uncertainties in the DIRBE absolute 
calibration. 

3.3 The Small Magellanic Cloud (SMC) 

How successful is this dust model in reproducing the emission observed from 
other galaxies? The SMC, with metaUicity -10% of solar (Dufour 1984), and 
dust-to-gas ratio ~10% of the Milky Way (Bouchet et al. 1985) is a good 
test case. Li & Draine (2002) sum over both quiescent medium and regions of 
active star formation in the SMC, find that to match the observed spectrum 
they require a distribution of starlight intensities. Following Dale et al. (2001) 
they adopt a power-law distribution of starlight intensities. 
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Fig. 22. IR surface brightness toward (a) I = 44°20', b = -0°20'; (b) / = 47° 45'. 
Diamonds: DIRBE photometry. Thin line, upper panel: 5-12 /im spectrum measured 
by MIRS on IRTS (Onaka et al. 1996). Thin line, lower panel: 2.8-3.9 fim spectrum 
measured by NIRS on IRTS (Tanaka et al. 1996). Solid line: carbonaceous-silicate 
dust model illuminated by 2xISRF, plus reddened starlight. Prom Li & Draine 
(2001). 



The dust extinction in the SMC is known to differ from Milky Way dust, 
primarily in the absence of a 2 175 A extinction bump on most (but not all) 
SMC sightlines (Gordon & Clayton 1998). Since we take the view that the 
2 175 A feature is due to PAHs, the absence of this feature implies a low PAH 
abundance. We adopt the dust mixture of Wcingartner & Draine (2001a), 
which reproduces the measured SMC bar extinction curve. 

With suitable choice of illuminating radiation field, the dust model is able 
to reproduce the observed IR emission from the SMC. It appears that the 
carbonaceous/silicate grain model (Weingartner & Draine 2001a; Li & Draine 
2001) can reproduce both the observed interstellar extinction by dust and the 
observed IR/FIR emission from dust. 




Fig. 23. Comparison between dust model (solid line) and the observed spectrum of 

the SMC obtained by COBE/DIRBE (diamonds) and IRAS (squares) averaged over 
a 6.25 deg^ region including the optical bar and the Eastern Wing (Stanimirovic et 
al 2000). Model spectrum includes infrared emission from dust heated by a range 
of starlight intensities, plus starlight, which dominates at A ^; 7 fim. Triangles show 
the model spectrum convolved with the DIRBE filters. Taken from Li & Draine 
(2002). 



4 Charging of Interstellar Dust 

Why do we care about the charging of dust grains? There are a number of 

reasons: 

• Charged dust grains are coupled to the magnetic field, whereas neutral 
grains are not. This is important not only for the motions of the dust 
grains, but also as a mechanism for couphng magnetic fields to neutral 
gas. 

• Charged dust grains undergo stronger gas drag due to Coulomb interac- 
tion with ions in the gas. 

• The photoelectric charging process injects energetic photoelectrons into 
the gas, which is a major mechanism for heating interstellar gas. 

• Neutral and negatively-charged dust grains can play an important role 
in neutralization of ions. 
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• In cold, dense regions, the dust grains may "lock up" a substantial frac- 
tion of the electrons. 



4.1 Collisions with Electrons and Ions 

Consider a spherical grain of radius a and charge ZgC, and a charged particle 
of charge ze, kinetic energy E, on a trajectory with impact parameter b. If the 
interaction potential is takcm to be a Coulomb potential, then conservation of 
energy and angular momentum allow one to find the critical impact parameter 
6crii for the particle to just graze the grain surface. Trajectories with b < 6crit 
will collide with the grain, and trajectories with b > 6crit will miss. The 
collision cross section is simply 



C = 7r6„it = na 



1 _ 



for £; > zeU , (59) 
= for E< zeU , (60) 

where U = Zgc/a is the electrostatic potential at the grain surface. If the 
velocity distribution is thermal, then the rate at which the charged particles 
reach the grain surface is just 



1/2 

= n.irn^ I I 

dt 



dN ^(ikrV'^ r^^zeu/kT iorzeU>0 



This is a wonderfully simple result. If the sphere is located in a plasma of 
electrons and ions, and the only charging process is coUisional, then the sphere 
will become negatively charged since the electrons move more rapidly than 
the ions. The steady-state charge is such that suppression of the electron 
arrival rate, and enhancement of the ion arrival rate, make them equal: 





r 




\-Kmi ^ 







(62) 



where we have assumed Zj = +1. Taking logarithms of both sides we obtain 



f^ = -llnf!^^UlnfflVl-fl--) 
kT 2 \mej \sej \ kTj 



(63) 



with the solution eU/kT = —2.51 for a hydrogen plasma {mi/nie = 1836) 
with Si = Se- 

However, the Coulomb potential does not fully describe the interaction be- 
tween a charged particle and a grain. In the case of a neutral grain, the electric 
field of the approaching charged particle will polarize the grain, resulting in an 
attractive potential - the same effect that causes ion-neutral scattering rates 
to be large, even at low temperatures. ^ Draine & Sutin (1987) have calculated 

^ This is just the familiar "image potential" from electrostatics for a charge near 
a conductor. 
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collision rates including the image potential. Image potential effects are im- 
portant for neutral grains when e^/a kT, or aT ;$ 1.67x 10^^ cmK. Thus it 
is an important correction for small grains in cold gas (e.g., a ^ 1.67x 10^*^ cm 
for T = 100 K). 

The "sticking coefficients" for small electrons and ions are not well known. 
In the case of an ion with an ionization potential larger than the work function 
for the grain material, it seems likely that the ion, upon arrival at the grain 
surface, will either remain stuck to the grain or will strip an electron from the 
grain and depart; in either case the ion "sticking coefficient" s/ = 1. The case 
of impinging electrons is less clear. The electron might be elastically reflected 
from the grain surface, or the electron might pass through a small grain and 
out the other side without energy loss, in which case the electron would have 
enough energy to return to infinity. 

Experimental sticking coefflcients for small neutral molecules or molecular 
ions can be estimated by dividing the measured electron capture rate coef- 
ficient by the estimated rate at which the electrons would reach the surface 
of the ion. Figs. [21 and [23 show the electron sticking coefficients for small 
neutral and charged carbonaceous clusters containing 10 < Nc < 100 carbon 
atoms. 

In hot gas, it is possible for impinging electrons to eject an electron from 
the grain. For molecules, this is called "ionization" ; for larger objects, it is 
called "secondary electron emission". It is therefore possible for grains to 
become positively charged by collisions in a hot (T ^ 10^ K) plasma (Draine 
& Salpeter 1979a). 

4.2 Photoelectric Emission 

The electrons in a neutral grain can be thought of a being confined witliin 
a potential well, as shown in Fig. I^Hl The potential well is produced by a 
"double layer" of charge at the boundary of the grain (the electron charge 
density extends beyond the ion charge density, resulting in a region near the 
surface with a radial electric field - see standard texts on solid state physics, 
e.g., Ashcroft & Mermin 1976). 
Three cases are shown: 

• For a grain with Z > 0, the excited electron will be subject to a long-range 
Coulomb attraction by the rest of the charge. In this case the electron 
which physically gets outside of the grain must still have additional energy 
in order to reach oo. 

• For a grain with charge Z ~ —1, the excited electron sees a system with 
zero net charge, so once the electron is physically outside the double layer, 
it will continue to infinity. 

• For a grain with charge Z < —1, an electron which gets outside the 
double layer will be accelerated away from the grain. An excited electron 
can in principle penetrate the double layer by tunneling, so that it is not 
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Fig. 24. Electron sticking coefficient for small neutral grains. From Weingartner & 
Draine (2001c). 

actually necessary for the electron to have an energy higher than the top 
of the double layer. 

An insulating grain has all the "valence band" energy levels occupied when it 
is neutral. The difference in energy between the top of the valence band and 
the first available vacant energy level (in the "conduction band" ) is referred 
to as the "band gap". If an insulating grain is to be negatively charged, 
the excess electrons must occupy available energy levels in the conduction 
band. Less energy is required to photoeject one of these "excess" electrons 
(a process known as "photodetachment" ) than to eject an electron from the 
top of the valence band. 

When an ultraviolet photon is absorbed in a dust grain, it raises an elec- 
tron to an excited state. If the electron has an energy high enough to reach 
infinity, and does not lose this energy to inelastic scattering, it will escape 
and be counted as a "photoelectron" . To calculate the rate of photoelectric 
charging of a grain, we require the "photoelectric yield" Y{hv, Z, a), which is 
the probability that absorption of a photon hv will produce a photoelectron. 
To calculate the rate of photoelectric heating of the gas, we also need to 
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Fig. 25. Electron sticking coefficient for small grains with charge Z = +1. From 
Weingartner & Draine (2001c). 

know the wicrgy distribution function f{E, hv, Z, a) of the energy E of the 
photoelectrons leaving the grain. 

The photoelectric ejection process consists of four stages: 

1. Excitation of an electron of appropriate energy. 

2. Motion of the electron from the point of excitation to the grain surface, 

3. Penetration of the surface layer and overcoming of the image potential. 

4. Once outside the grain surface, the electron must overcome the Coulomb 
potential (if the grain is now positively charged) to reach infinity. 

Weingartner & Draine (2001c) have estimated the photoelectric yield for 
small grains, both neutral and charged, taking into account the above effects. 
Allowance is made for dependence of the image potential on the grain radius 
a, and the possibility of tunneling through the double layer in the case of a 
negatively charged grain. 

Watson (1973) pointed out that very small particles should have enhanced 
photoelectric yields because of step 2 above: in a small particle, photoexcita- 
tions will on average be closer to the grain surface, increasing the probability 
that the excited electron will reach the grain surface without energy loss. 
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Fig. 26. Model for the potential confining electrons in a grain with charge Ze. 
Shaded regions show occupied energy levels. W is the work function for bulk ma- 
terial; IP{Z) is the ionization potential; EA(Z) is the electron affinity; hvpct is the 
photoelectric threshold energy; hupdt is the photodetachment threshold energy (the 
minimum photon energy required to detach an electron from a negatively charged 
grain); IPv{Z) is the energy difference between infinity and the top of the valence 
band; -Ebg is the bandgap energy for the material (-Ebg = for a metal). Taken 
from Weingartner & Draine (2001c). 



Draine (1978) proposed a simple function to estimate the magnitude of this 
yield enhancement, which we use: 



a 



-2a + 2- 2exp(-a) 
/32-2/J-|-2-2exp(-/J) 



a 



a = — + — 



(64) 
(65) 



where la = 
and L, the 



A/[47rIm(m)] is the photon attenuation length in the material, 
"electron escape length" , is essentially the mean free path against 
energy loss to inelastic scattering for an excited electron. In the limit a <C 1 
we have 

a la + le 

= ■ 



(66) 
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hi^(eV) 

Fig. 27. Photoelectric yield for neutral coronene, as measured by Vcrstractc et 
al. (1990) and as estimated in the approach of Weingartner & Draine (2001c). Eq. 
11, 12, 13, 16 refer to equations in Weingartner & Draine (2001c), from which this 
figure is taken. 

Photon attenuation lengths la ~ 300A are typical in the vacuum ultravi- 
olet. Martin ct al. (1987) report le ~ 9 A for 6 eV electrons in thin carbon 
films, and McFeely et al. (1990) report le ~ 6A for 8 eV electrons in Si02. 
Weingartner & Draine (2001c) adopt Ig = lOA, independent of energy, for 
both graphite and silicate grains. Since 1^ <C la, it is clear that very large 
yield enhancement factors yi are possible for very small grains. 

Weingartner & Draine write 

Y{hi/, Z, a) = min [yo x yi(hv, a)] x y2(hv, Z, a) (67) 

where 2/2 is the fraction of electrons which, having crossed the grain surface, 
have sufficient energy to overcome the long-range coulomb attraction if Z > 
(for Z < 0, y2 = 1). For PAHs and larger carbonaceous grains, yo is chosen 
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Fig. 28. Photoelectric yields (photoelectrons per absorbed photon) for neutral 
graphite and silicate grains as a function of incident photon energy hi^, for selected 
values of the grain radius a. From Weingartner & Draine (2001c). 

SO that Y = ?/o2/i2/2 approximates the photoionization yield measured for 
coronene. 

Fig. 1281 shows the photoelectric yields estimated for neutral carbonaceous 
and silicate grains of various radii. 

4.3 Charge Distribution Functions 

Fig. 1291 shows grain charge distribution functions calculated for dust in a 
molecular cloud region. The H nucleon density is taken to be tt-h = 300 cm^"^, 
the gas temperature T = 25K, and the fractional ionization ng/nn = 10^**. 
The radiation field is assumed to have the spectrum of the MMP ISRF, but 
with intensity a factor of 10 below the MMP value. Under these conditions 
small grains (a ^ 200A) remain negatively charged. For larger grains, the 
reduced importance of "image charge" effects (as well as a slight reduction 
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Fig. 29. Charge distribution functions for carbonaceous grains and silicate grains 
in cold molecular gas. 
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Fig. 30. Charge distribution functions for carbonaceous grains and silicate grains 
in "cold neutral medium" diffuse HI. 
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in the work function) tilts the balance in favor of photoelectric charging, but 
the grain potentials remain very small. 

Fig. 1301 shows grain charge distribution functions calculated for a dif- 
fuse region with starlight equal to the MMP ISRF. The electron density 
rie = 0.03 cm~'^ is taken to be the same as for the molecular cloud region 
of Fig. but the increased starlight intensity allows photoelectric emis- 
sion to dominate the charging for grains with a > 25A. function The "CNM" 
conditions are appropriate to diffuse HI at 100 K. Despite the presence of ul- 
traviolet radiation producing photoelectric emission, note that an appreciable 
fraction of the small grains can be neutral or even negatively charged. 

We see that the competition between collisional charging and photoelec- 
tric charging can go either way, depending on the grain composition, the 
grain size a, the electron density rig and temperature T, and the spectrum 
and intensity of the ultraviolet background due to starlight. 

5 Dynamics of Interstellar Dust 

What are the velocities of interstellar grains? To answer this question, we 
must understand the forces acting to accelerate or decelerate the grains. 



5.1 Gas Drag 



If the grain is moving relative to the gas, it will experience a drag force 
as momentum is transferred from the grain to the gas. This drag force is 
approximately -Fkrag = MgrVgr/rgr where Mgr is the grain mass, Vgr is the 
grain speed relative to the gas, and Tgr is the time for the grain to collide 
with its own mass in gas particles. For the case of an uncharged spherical 
grain in neutral gas, the problem can be solved exactly, and approximated to 
within 1% by 



Qtt 
64' 



1/2 



rriiV 
2kT 



(68) 



where the sum is over the species (H, H2, He) in the gas (Draine & Salpeter 
1979a). If the grain is charged and ions are present, there is an additional 
"Coulomb drag" term (Draine & Salpeter 1979a), but this is numerically 
unimportant in regions of fractional ionization ^ 10~^. Eq. (|68|l is obtained 
for either perfectly inelastic collision (where the gas atom is assumed to 
"stick") to the grain or perfectly elastic specular reflections. For subsonic 
motion, the drag time is 



MgrV pa 



rag 



(69) 
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For a cool region where the hydrogen is primarily molecular, 
Tdrag = 1.2 X lO^yr 



300cm-3\ /25K^ 



^0.1 fj,m. J \3gcm J \ J \ T 

while in a more diffuse region where the hydrogen is atomic, 



(70) 



^'^-^ = 'ioAj;^) (s^j [^) • ^''^ 

5.2 Lorentz Force 

Dust grains are charged, and there are magnetic fields in clouds. The cy- 
clotron period is 

2n ^ 2TrMgrC 

UJc \{Z)\eB ^' ' 

5.2xl0''yr/ a \3 / p \f3fiG\ 



\{Z)\ Vl0-5cm/ V3gcm-3 



where we take (Z) because the grain charge fluctuates on a time short com- 
pared to the cyclotron period. Values of {Z) are indicated in Figs. [5^1 and 1201 
The character of the grain dynamics depends on the ratio of the drag force 
to the Lorentz force, or 



3 \{Z)\eB TMfcT 
- 53 m\ (^) ri^V f^V' C-^) (75) 



3/iG /\ a J \ T J \?iH 

where we have assumed the hydrogen to be atomic. For grains in cold clouds, 
we see that we usually have cOcTdrag 3> 1: grains are strongly coupled to the 
magnetic field, except for grain sizes where collisional charging and photo- 
electric charging happen to just balance so that \{Z)\ ^ 0.02(a/10^^ cm)^. 



5.3 Radiation Pressure 



The interstellar radiation field is in general anisotropic. As a result, dust 
grains are subject to "radiation pressure" forces which can cause drift of 
the grain relative to the gas. This drift can be important as a means of 
transporting dust through gas, but also because the drifting grains transfer 
momentum to the magnetic field (via the Lorentz force) and to the gas (via 
gas drag). 

Consider a unidirectional beam of photons of energy hiy, with energy 
density Ui^dv in interval {v, v -\- dv). The photons have momentum hv /c^ and 
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this momentum can be transferred to the grain by absorption or by scattering. 
It is easy to see that the rate of momentum transfer is 

F,ad - / [Cabs + Csca(l " (cOS 9))] — (76) 

J hu c 

^ TTO^ j U^dv Qp^{v) (3pr(i^) = Qabs + Qsca(l - (C0S6')) . (77) 

where Qabs = Cabs/^ra^, Qsca = Csca./'^o,'^ , and (cos0) is the mean of the 
cosine of the scattering angle. 

We are generally interested in values of Qabs and Qpr averaged over the 
spectrum of the radiation. 

We are often content to approximate the radiation from a star by a 
blackbody spectrum. In neutral regions of the interstellar medium we are 
frequently interested in averages over a blackbody spectrum with a cutoff at 
13.6 eV, and it is useful to calculate spectrum-averaged values 

_ j;-d. B^{T)Q{.) 
Jo diy B^{T) 

where B^{T) is the Planck function H53|l and /ifmax — 13.6 eV. 

We will also be interested in Q values averaged over the interstellar ra- 
diation field, for which we use the estimate of Mathis, Mezger, & Panagia 
(1983). 

Suppose that at some point in the interstellar medium the starlight has a 
net energy flux cZ\itrad- By summing over the Skymap Star Catalog (Slater & 
Hashmall 1992), the Tycho Catalog (ESA 1997), and the all-sky ultraviolet 
observations of the S2/68 experiment on the TD-IA satellite, Weingartner & 
Draine (2001b) estimated the anisotropy in the starlight background at V, 
B, and 4 ultraviolet bands. In Fig.|221we show the net dipole moment of the 
radiation field at these 6 wavelengths. 



5.4 Recoil from Photoelectric Emission and Photodesorption 

Radiation pressure is, of course, just the transfer of photon momentum to 
the grain. However, radiation can serve to energize other processes which can 
exert a thrust on the grain: photoelectric emission and photodesorption. To 
see that these might be important, it is sufficient to compare the momentum 
of an hy « 10 eV photon with the momentum of the photoelectron it might 
produce, if the photoelectron kinetic energy = / x hv: 

p(electron) _ {2m,fhvf'^ -ZIQJ^'^ flOcYV'^ ^ ^^^^ 
p(photon) hv I c \ hv J 
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Fig. 31. Radiation pressure efHciency factor (Qpr) for neutral carbonaceous and 
silicate grains, averaged over the interstellar radiation field (ISRF) and blackbody 
spectra with indicated color tempertures, cut off at 13.6 eV. Taken from Weingart- 
ner & Draine (2001b). 



Similarly, suppose that a photon caused an H2 molecule to be desorbed with 
kinetic energy f x hv: 



p(photon) hv/c \ hv 

Therefore if an appreciable fraction of absorbed photons produce photoelec- 
trons or photodesorb atoms or molecules, there could be a significant recoil 
on the grain. 

For a spherical grain in an isotropic radiation field, the photoelectrons 
would be emitted isotropically, with zero net thrust. Suppose, however, that 
we have a unidirectional radiation field (e.g., light from a single star). In this 
case we might expect a higher rate of photoelectron emission from the "bright 
side" than from the "dark side" of the grain. However, very small grains do not 
have a "dark side" - they are effectively transparent to the incident radiation, 
so one must try to estimate how this depends on grain size. How can we 
estimate the bright side/dark side asymmetry? A simple model is to assume 



Astrophysics of Dust in Cold Clouds 



53 




Fig. 32. Dipole fraction for the interstellar radiation field. Upper curve is the ratio 
of the measured dipole moment to the measured starlight (the measurements are 
incomplete). Lower curve is the ratio of the measured dipole moment to the total 
estimated ISRF. The true fraction is probably between these two curves. After 
Weingartner & Draine (2001b). 



that the rate of photoelectric emission from a point on the grain surface is 
proportional to |Ep just inside the surface. Using Mie theory to calculate 
the electric field intensity E in a spherical grain, one can then evaluate the 
relative rates of photoelectric emission over the grain surface. Using this 
approach, Weingartner & Draine (2001b) have investigated the net thrust due 
to escaping photoelectrons - equal and opposite to the net momentum which 
the escaping electrons have when they reach infinity. The photoelectric yield, 
and the energy of the photoelectrons, is of course affected by the charge on the 
grain, and therefore depends on the environmental conditions. The analysis 
is complicated by the fact that if the grain is charged, photoelectrons with 
nonzero angular momentum will escape to infinity in directions different from 
the direction they are travelling when they emerge from the grain surface 
(until it reaches infinity, the electric field of the photoelectron exerts a force 
on the grain). 

Fig. 1331 shows the results of Weingartner & Draine (2001b) for the ratio of 
the photoelectric thrust to the radiation pressure force, as a function of grain 
size. We see that this is an order-unity correction for grains with a « lOOA 
in the CNM. 

For photoelectric emission, it is reasonable to assume that the grain is 
"spherically symmetric" prior to arrival of the photon - there are just as many 
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a (A) 

Fig. 33. Photoelectric tlirust and pliotodesorption tiirust relative to radiation pres- 
sure for carbonaceous and silicate grains as a function of grain radius a. Curve 
labelled "CNM" is for cold diffuse clouds with = .03 cm"^, T = 100 K, and the 
standard ISRF. Note that the photodesorption force can be negative for grains with 
radii lOOA <a< 500A. From Weingartner & Draine (2001b). 

electrons available for photoejection on the bright side as on the dark side of 
the grain. For photodesorption, however, this is not a reasonable assumption. 
As we will see in below, grains are expected to be spinning rapidly, with 
their principal axis of largest moment of inertia expected to be aligned with 
the angular momentum. Because of this gyroscopic stabilization, one may 
expect that if photodesorption is a rapid process, then the "bright" side of 
the grain will have fewer adsorbed molecules ready to be photodesorbed, an 
effect which will depend on the orientation of the grain's spin axis with the 
direction of radiation anisotropy. 

If photodesorption is rapid, the thrust will be limited by the need to 
resupply the grain surface with new molecules to be photodesorbed. The 
story has a number of additional complications: 

• The result will depend on the wavelength-dependence of the photodesorp- 
tion cross sections for the adsorbed species. This is essentially unknown. 
Weingartner & Draine simply assume that the photodesorption rate for 
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a molecule on the grain surface is proportional to |Ep just outside the 
grain surface, integrated from 6-13.6 eV. 

• Adsorbed species may diffuse over the grain surface, so that the "bright 
side" could be replenished from the "dark side" . 

• If photodesorption is rapid, the surface coverage of adsorbates will de- 
pend on the orientation of the grain spin axis relative to the direction of 
starlight anisotropy. 

• The spinning grain can undergo "polar flips" - the "thermal flipping" 
process described by Lazarian & Draine (1999a). 

• If the grain is moving, the grain motion will alter the rates at which 
molecules are resupplied to the different parts of the grain surface: the 
"bright side" may suffer a reduction in the rate of arrival of atoms or 
molecules from the gas. 

• H atoms which are not photodesorbed might react on the grain surface 
to form H2, and could recoil from the grain with substantial kinetic en- 
ergy. This process might happen more rapidly on the dark side where 
the H concentration is higher, thus acting to oppose the thrust due to 
photodesorption. 

With so many uncertainties, it is not possible to reach definite quantitative 
estimates for the photodesorption force, but Weingartner & Draine (2001b) 
evaluate one plausible set of assumptions; results are shown in Fig. 123 for 
grains in the "warm neutral medium" and the "cold neutral medium" . It is 
curious that there is a range of grain sizes (100 - 400A for carbonaceous grains, 
100 - lOOOA for silicate grains) where they find that the photodesorption force 
is negative - this occurs because for this range of sizes, interference effects 
cause |Ep to be larger on what would have been expected to be the dark side 
of the grain if the grain were opaque and geometric optics were applicable. 
Note that the larger grains do behave as "expected" for macroscopic targets. 

The "bottom line" is that photodesorption can make a significant contri- 
bution to the net force on the grain, which (coincidentally) is expected to be 
comparable in magnitude to the radiation pressure force and the photoelec- 
tric thrust, but which can in principle be of the opposite sign, depending on 
detailed physics of photodesorption which are not known at this time. 

5.5 Drift Velocities for Interstellar Grains 

For the special case where the radiation anisotropy is parallel to the magnetic 
field direction, we can obtain the steady-state grain drift velocity by balancing 
the drag force against the radiation-related forces: 

-Fdrag — -Fr-ad + ^pc + ^pd , (82) 

where Fpc and Fpd are the forces due to photoelectric emission and photodes- 
orption. If we assume subsonic motion, and neglect Coulomb drag, this can 
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Fig. 34. Drift velocities for carbonaceous and silicate grains, as a function of grain 
size, in the CNM and WNM, for a radiation anisotropy zAwrad/uisRF = 0.1. Taken 
from Weingartner & Draine (2001b). The minimum in drift velocity for a ~ 250A 
carbonaceous grains or a « 350A silicate grains in the CNM is due to the reversal 
of the photodesorption force for grains in this size range (see Fig.Onj- 
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be written 

^^drift — f ^^fatl ^ "'^ ^ Fpc/Frad + -Ppd/ -Frad ^gg^ 

so we see that the grain will drift with a "Mach number" which is approxi- 
mately equal to (Qpr) times the ratio of the anisotropic radiation pressure to 
the gas pressure. Since Qpr reaches values of order unity (see Fig. 131(1. we may 
expect peak Mach numbers or order 0.2 or so - the ratio of the anisotropic 
radiation energy density to the gas pressure in the diffuse ISM. 

In Figure El we show drift speeds calculated by Weingartner & Draine 
for carbonaceous and silicate grains in the WNM and CNM. The drift speeds 
are clearly a strong function of grain size, and appear to be negligible for 
small a < lOOA grains. For larger grains the drift velocities are modest but 
potentially significant since they could be sustained for long times: in the 
CNM, a radiation anisotropy of only Z\urad/urad = 0.1 would give a 0.3 /im 
carbonaceous grain a drift speed of ~ 0.04kms~^, so in 10 Myr it could drift 
0.4 pc, and in the WNM the drift speeds are about an order of magnitude 
larger. This could in principle lead to removal of the large grains from some 
gas elements and concentration in others. 



6 Rotational Dynamics of Interstellar Dust 

The rotational dynamics of interstellar grains constitutes a fascinating story, 
not yet fully told. It is intimately tied to the long-standing problem of inter- 
stellar grain alignment, which we will discuss in |7| 



6.1 Brownian Rotation? 

Early discussions of the rotational dynamics of interstellar grains generally 
assumed the grains to be rigid spheres, which appeared to capture the essen- 
tial physics, and discussed the scattering of impinging atoms from the grain 
surface. As far as grain rotation was concerned, it seemed obvious that the 
grain would simply undergo a form of Brownian motion, exchanging linear 
and angular momentum with the gas. Once statistical equilibrium was at- 
tained, a rigid grain (with purely elastic scattering of atoms from its surface) 
would have an expectation value of fcTgas/2 for the kinetic energy in each 
of its degrees of freedom: three translational and three rotational. Thus the 
r.m.s. translational velocity of the grain would be 



1/2 / ni.rr \ 1/2 



2,1/2 _ fikTV'^ _ f 9kT 



s VlOOK 
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which is negligible under all foreseeable circumstances. The r.m.s. rotation 
rate 



(86) 



4 , T \^/V3gcm-3V/VlO"'cm\'/' , 
4-6x104Hz — — — (87) 



100 K 

is more impressive, though it simply corresponds to an equatorial velocity 
equal to the r.m.s. translational velocity times a factor a/5/2 — 1.58. 



6.2 Suprathermal Rotation 

Following the discovery of grain alignment in 1949, theoretical studies of the 
rotational dynamics were carried out to try to understand the alignment 
process, including the seminal paper by Davis & Greenstein (1951). It was 
recognized that, since interstellar grains were likely to have vibrational tem- 
peratures differing from the gas temperature, eq. H87|) needed to be modified, 
and it was understood that for a spherical grain the temperature T appearing 
m eq. l|57|l should be some weighted mean of the gas temperature and the 
grain temperature, since collisions of atoms with the grain surface will not 
be perfectly elastic. 

However, in 1975 - 24 years after Davis & Greenstein's study - Purcell 
(1975) realized that interstellar grains could act as "pinwheels". After all, a 
simple pinwheel weighing many grams can readily be put into rotation with 
a rotational kinetic energy many orders of magnitude greater than kT by the 
air molecules which excite it. Purcell (1979) showed, in fact, that there were 
3 distinct physical processes acting on interstellar grains which could each, 
acting alone, cause a grain to rotate with a rotational energy (l/2)/ci;^ 3> kT: 

1 . The effective "accomodation coefficient" will probably vary over the grain 
surface, due to chemical inhomogeneities and/or geometrical irregulari- 
ties. Such variations in accomodation coefficient on a nonspherical grain 
will lead to a systematic torque (i.e., one which whose time-average is 
not expected to be zero) if the gas temperature and grain temperature 
differ. 

2. Photoelectrons ejected from the grain may, on average, carry away angu- 
lar momentum, resulting in a systematic torque. 

3. H2 formation on the grain surface, with the nascent H2 molecules ejected 
with a significant kinetic energy, is expected to result in a large system- 
atic torque, particularly if the H2 formation occurs at a relatively small 
number of "active sites" on the grain surface. 

Purcell showed that of the three processes, the H2 formation process was 
likely to be the most important in regions where the H is atomic. 
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I6.2L 1 Suprathermal Rotation Driven by H2 Formation 

It is instructive to examine the argument for suprathermal rotation driven 
by H2 formation at active sites on the grain surface. 

To simplify the argument, let us consider a region where only atomic H is 
present (no He), and consider a cube of sides 26 x 26 x 26, with volume 86"^, 
surface area 246^, and effective radius acs — (6/77)^/^6 — 1.246. 

Choose a site at random on the grain surface, and assume that H2 molecules 
are formed at that site at a rate A^i . For simplicity, suppose that the newly- 
formed H2 molecules leave normal to the surface with kinetic energy En^ . 
The expectation value for the square of the moment arm is then just (2/3)6^, 
and the expectation value for the mean square torque exerted by this site is 

= N^[mii,vii,f X ^6^ 

= ^N^maEH.b^ . 

Now suppose that there are a total of NgUe such sites on the grain surface, 
and suppose that a fraction 7 of the H atoms which impinge on the surface 
are converted to H2 (so there are 7/2 H2 formation events per arriving H 
atom). In a thermal gas, the rate per area at which atoms collide with a 
convex surface is just n(kT/2Trmy^'^, so 



Since the sites are assumed to be randomly located, the torques add like a 
random walk, so that the expectation value for the square of the total torque 
is just 

(rl) = NsUeF^ (91) 

1 Q2 1 

= — ^^n{RfE^,kT b"^ . (92) 

The rotating grain will experience gas drag. For slow rotation rates u, it is 
not difficult to show that the drag torque is 



1/2 

^drag = -Y^iH) {^^^ ] ■ (93) 



80 f m-akT\ 



To find (w^), we now set {Fl^) = {Ej^^ ) to obta 



2, 27 Y En, 1 



50 Nsite TOH 6^ 



or a rotational kinetic energy, relative to the thermal value, 

/(^') ^ 3 M 7^ En, 
3kT 25 mn NgUe kT 



(95) 
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— "£:(T^)'(^)(!f)(T)-<-) 

where M = p is the grain mass. Now this is no small number! The nmiiber 
of active sites is micertain. If we suppose that there is a surface area per 
active site, then NgUe — and 

(97) 



- 8.4x lOV I r 1 I s 1 I ? I I 1 I 

The r.m.s. rotation rate is 



^-.(^)(l^)(g.) , .s, 

so that we might have MHz rotational frequencies for 6 < 3 x 10^^ cm! It is 
apparent that systematic torques can play a major role in grain dynamics. 

Each II2 formation site can be thought of as being Uke a small rocket 
thruster attached to the grain surface: the systematic torque due to II2 for- 
mation is fixed in body coordinates, so long as the H2 formation sites do not 
change. The kinetics of H2 formation on grain surfaces is poorly understood 
at this time, and it is not certain how long-lived the active sites are hkely to 
be. 

ini2l2 Radiative Torques Due to Starlight 

The three torques identified by PurceU do indeed appear to be individuaUy 
capable of driving grains to suprathermal rotation. In fact, there appears to 
be an additional physical process which can compete with the H2 formation 
torque for grains larger than ~ lO^^ cm: torques exerted on interstellar grains 
by starlight. 

We discussed above the "radiation pressure" force which anisotropic starlight 
could exert on a grain. It turns out that starlight can also exert a torque on 
grains, and these torques can be dynamicaUy important! 

How can we address this problem? Analytic progress is difficult, because 
the only cases for which we have analytic solutions are cases where radia- 
tive torques vanish: (1) grains of complex shape but small compared to the 
wavelength; (2) spheres; (3) spheroids. In the first case the torques vanish 
because only the dipole response of the grain is important, and the dipole 
radiation does not carry angular angular momentum. In the case of spheres 
and spheroids, there is no torque because of the symmetry of the target. 

In order to obtain a nonzero torque we must have an asymmetric target. 
Draine & Weingartner (1996) discussed use of the discrete dipole approxi- 
mation to calculate radiative torques. The DDSCAT code (Draine & Flatau 
2000) now provides the capabihty to compute torques. For an asymmetric 
target illuminated by a beam of unpolarized radiation, the direction and 
magnitude of the torque depends on the orientation of the target relative to 
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b) 



Fig. 35. Macroscopic examples of targets which would be subject to radiative 
torques. From Draine & Weingartner (1996). 



the beam. It is useful to define a "torque efRciency vector" Qr defined so 
that the torque on the grain is given by 



2 ^ ^ 



(99) 



where the effective radius Ooff is the radius of an equal vohimc sphere. We 
will generally want to average Qr over some spectrum; Draine & Weingartner 
(1996) calculate 



(Qr) 



J XuxdX 



- ^ / XuxdX 
J u\dX 



(100) 
(101) 



for MA taken to be the MMP ISRF (for which A = 1.202 /xm). With these 
definitions, the torque 



Trad = 7raeff"rad7^(Qr) 
ZTT 



(102) 



We recall from classical mechanics that an arbitrary object has three 
principal axes, ai, a2, §.3, and the moment of inertia tensor is diagonal in a 
coordinate system aligned with these three axes. Let the eigenvalues of the 
moment of inertia tensor be Ii > I2 ^ h- 
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For purposes of discussion, let us assume that Ii > I2 = h- K the grain 
has angular momentum J, then the rotational kinetic energy is 

(Jcos/3)^ , (Jsin/3)^ 

= 2/1 + 

= ^ + J^^^^Bin^/3 , (104) 

with the rotational kinetic energy minized for /3 = 0. As we will see below, 
when J is large we expect the grain to remain close to this minimum energy 
state, with sin^ /3 ^ 1. 

The grain spins rapidly, so we are concerned with the torque after aver- 
aging over rotations of the grain around the axis ai. This rotation-averaged 
torque therefore depends only on a single angle - the angle between ai 
and the direction of the incident radiation. 

We will see below that we expect a spinning dust grain to precess around 
the direction of the local magnetic field B, so what we really need is the torque 
averaged over grain rotations and over precession of the axis ai around the 
magnetic field direction. For a given grain and spectrum of radiation, this 
will now depend on two angles: the angle ip between B and the direction of 
the radiation, and the angle f between ai and B. 

For purposes of spinup (or spindown), we are interested in the projection 
of the torque efficiency vector along the spin axis, i.e., the axis of largest 
moment of inertia: Qr(0) • a-^. Draine & Weingartner (1996, 1997) define 
ff(^, (p) = (Qr ■ ai) averaged over rotation around ai and precession of ai 
around B. They calculate i?(C, 0) for three grain shapes with effective radii 
floff = 0.2 /im. The numerical values of H of course depend on ^ and 0, but a 
typical value of \H\ would be \H\ « 0.005 (see Fig. 5 of Draine & Weingartner 
1997), so that the spinup torque would be 

{Trad ■ ai) = HXAUradT^alg . (105) 

It is of interest to compare this to the spinup torque due to H2 formation. 
We can rewrite eq. replacing b by {n/dY^^aeS, and project onto a random 
direction ai to obtain 

Thus the ratio of the radiative spinup torque (|105|l to the r.m.s. value of the 
H2 formation torque is 



{Trad ■ ai) HX AUrad 



((Fe, • ai)2)i/2 O.Wl^iEuJkTy/^s n(R)kT 



2.92 ( H \ (J_\'^ (^\'^ (]^\ ^""r-d 
7 V0.005y VIOOkJ \EnJ \ s J n{ll)kT 



(107) 
(108) 
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Adopting the 10% anisotropy estimated by Weingartner & Draine (2001b), we 
have Aurad ~ 8.64 x 10"" crgcm-^. With n(H) = SGcni-^, and T = 100 K, 
we have Aurad/n^H^kT = 0.21. Therefore with an H2 formation efficiency 
7 < 0.5 the radiative torque would exceed the H2 formation torque if the 
other parameters have nominal vahics (in particular, s = lOA). Many of the 
numbers are uncertain; the important point is that for a^s = 0.2 /xm grains 
the radiative torque has a magnitude which is comparable to our estimate 
for the H2 formation torque. 

7 Alignment of Interstellar Dust 




Fig. 36. Linear polarization of light from 1800 stars. From Mathewson & Ford 
1970. 

Polarization of starlight was discovcTC'd in 1949. It was immediately de- 
duced that the only plausible mechanism was "linear dichroism" of the inter- 
stellar medium due to aligned dust grains, but understanding the mechanism 
responsible for this alignment has frustrated astrophysicists for half a century. 
In fact, the story of the efforts of astrophysicists to understand grain align- 
ment is humbling. It is not just because it has taken us so long to solve the 
problem the humbling lesson is that several very important and relatively 
simple physical effects were overlooked for decades. It was not that these 
physical effects were difficult to understand, or that they involved physics 
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which has only recently been understood. On the contrary - the physics in- 
volved would have been familiar to well-trained physicist in 1950. It was 
simply that important processes were overlooked: 

• The Rowland effect: a charged, spinning dust grain will develop a mag- 
netic moment due to its circulating charge (Martin 1971). 

• The Barnett effect: a spinning dust grain with unpaired electron spins 
will spontaneously magnetize (Dolginov & Mytrophanov 1976). 

• Suprathermal rotation due to dust-gas temperature differences (Purcell 
1975, 1979). 

• Suprathermal rotation due to photoelectric emission (Purcell 1975, 1979). 

• Suprathermal rotation due to H2 formation (Purcell 1975, 1979). 

• Viscoelastic dissipation of rotational kinetic energy due to time-varying 
stresses in a grain which is not rotating around a principal axis (Purcell 
1979) 

• "Barnett dissipation" of rotational kinetic energy due to the electron spin 
system (Purcell 1979). 

• Dissipation of rotational kinetic energy due to the nuclear spin system 
(Purcell 1979). 

• Suprathermal rotation due to starlight torques (Draine & Weingartner 
1996). 

• Fluctuation phenomena associated with Barnett dissipation (Lazarian & 
Roberge 1997; Lazarian & Draine 1999a) and coupling to the nuclear 
spins (Lazarian & Draine 1999b) 

Have we thought of all the important processes? Or are there other important 
phenomena to which our eyes are still closed? 

7.1 The Motion of Spinning Grains 

As seen above, dust grains spin very rapidly. Consider for the moment one 
spinning grain. If the angular momentum J does not happen to be along one 
of the grain's principal axes, then the grain's angular velocity oj will not be 
parallel to J. and will nutate around J. with a nutation rate that will be 
smaller than (but comparable to) the rotation rate. 

Now imagine observing a single grain as it spins at kHz rotation frequen- 
cies. Because the rotation and nutation are both rapid, we will quickly obtain 
a rotation-nutation average of the grain orientation. During this observation, 
the angular momentum J is conserved, and defines a preferred direction. 
Now consider light propagating in some particular direction. Because of the 
rotation-nutation averaging, the grain will tend to polarize the light with 
electric vector either exactly perpendicular to J, or exactly parallel to J, de- 
pending on whether the grain on average has its "long" axis perpendicular 
to J, or parallel to J. 
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This is for one grain. If a population of grains shows a tendency to po- 
larize light, then the individual grain angular momenta cannot be randomly- 
oriented - there must be a tendency for the angular momenta J to be aligned 
with some preferred direction in space. 

The problem of grain alignment is to account for the alignment of grain 
angular momenta J. 

Many alignment mechanisms have been proposed over the years. Some 
directly involve the galactic magnetic field (Davis & Greenstein 1951). Other 
proposed mechanisms rely on so-called "mechanical alignment" due to gas 
(Gold 1952) or cosmic rays (Salpeter & Wickramasinghe 1969) streaming 
relative to the grains. Harwit (1970) suggested that grain angular momenta 
might be aligned as a result of the angular momentum absorbed from starlight 
photons, if the starlight is anisotropic. It was only in 1971 that Martin (1971) 
pointed out that even in the case of nonmagnetic alignment processes, the 
observed alignment would be either exactly parallel or exactly perpendicular 
to the magnetic field, simply as a result of the magnetic moment on a spinning 
grain. 



7.2 Grain Magnetic Moments and Precession Around B 

Grain in diffuse regions tend to be electrically charged (see Fig. QUI) . If 
charged, a spinning dust grain will develop a magnetic moment due to the 
"Rowland effect" : the rotating charge on the grain constitutes a current circu- 
lating around a loop, and there will be a resulting magnetic field and magnetic 
moment interior to the loop. Recall that a current / flowing in a loop with 
area A generates a magnetic moment m = lA/c. The mean of x'^ + over 
the surface of a sphere is just (2/3)a^, so the magnetic moment due to charge 
Ze distributed uniformly over the surface of a spinning sphere is 

LuZe 2a^ uZea^ 

If u) is at an angle > Q relative to B, the spinning sphere will precess with 
a precession frequency 



jiB sin 9 5 ZeB 

LJp = 



Ilu sin 6 Stt pa^c 



(110) 



— (^) (^) ■<-) 

We have seen above (Fig. 1301) that in a CNM environment, we might expect 
an a = 10~^ cm grain to have (Z) w 10, for a precession period ~ 6000 yr. 
This time is short compared to the other dynamical times: the gas drag time 
H71|) . or (as we shall see below) the time for grain alignment to occur. 

We can therefore conclude that regardless of what process is responsi- 
ble for alignment of interstellar grains, each grain will have a time-averaged 
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angular momentum distribution which must be symmetric around B. There- 
fore the interstellar dust mixture will tend to polarize starlight either exactly 
parallel to B or exactly perpendicular to B, for any conceivable alignment 
mechanism! 



7.3 Barnett Effect 

Dolginov & Mytrophanov (1976) pointed out that the "Barnett effect" would 
result in grain magnetic moments much larger than due to the Rowland effect. 
The Barnett effect can be stated very simply: an uncharged object rotating 
with angular velocity lo tends to spontaneously magnetize, with a magneti- 
zation M = X'^/li where x is the susceptibility, and 7 is the "gyromagnetic 
ratio" for the material - the ratio of (magnetic moment) /(angular momen- 
tum) for the electron orbitals or electron spins responsible for the magneti- 
zability of the material. For electron spins, 7 = /Ub/(S./2) = e/2meC, where 
fiB = efi/2meC is the Bohr magneton. 

The spontaneous magnetization is equal to that which would be produced 
in a stationary sample in a (fictitous) applied magnetic field Hbc = w/7 = 
2meCU}/e - the so-called "Barnett equivalent field". 

At first sight the Barnett effect seems surprising, but it is not difficult to 
understand. Consider an unmagnetized sample (equal numbers of "spin-up" 
and "spin-down" electrons), with angular momentum J, spinning at angular 
velocity iv. At constant total angular momentum, the angular velocity of the 
sample can be reduced, with consequent reduction of the rotational kinetic 
energy if some of the angular momentum is put into the electron spin 

system. The electrons have magnetic moments, so with unequal numbers of 
electrons in the two spin states, the sample is now magnetized. The reduction 
in rotational kinetic energy is accompanied by heating of the grain, with an 
increase in entropy. The effect was observed in the laboratory by Barnett 
(1915). 

For normal paramagnetic materials, the magnetic susceptibility is just 

X = '>^■s^J''B/kTg, where is the number density of unpaired electrons, and Tg 
is the grain temperature. The magnetic moment due to the Barnett effect is 
thus 

^na^ nsiJ.% 2meC 

= ^ ■ (112) 

^ 3 kT e ^ ' 

The ratio of the "Barnett effect moment" to the "Rowland effect moment" is 
/x(Barnett) 2'Kansh'^ 



/i(Rowland) ZhTgirie 

2.8 X lOV « \ A20K 



Z V 10-5 cm/ \ Tg J VlO 



(113) 
(114) 



An unpaired electron density w 10^^ cm ^ is typical for normal para- 
magnetic materials, so it is clear that the Barnett effect is very important 



Astrophysics of Dust in Cold Clouds 



67 



for interstellar grains, and will result in precession periods of order days or 
months rather than the thousands of years for the Rowland effect. 

7.4 Alignment by Paramagnetic Dissipation: Davis-Greenstein 
MechEinism 

Davis and Greenstein (1951) described a mechanism which could produce 
grain alignment. Consider a grain spinning with angular velocity w in a static 
magnetic field Bq. Let 9 be the angle between Bq and lo. 

The grain material has a magnetic susceptibility Xi and it will try to 
magnetize in response to an applied magnetic field. In grain body coordinates, 
the applied field has a static component parallel to the grain rotation axis, 
with magnitude B|| = Bq cos 6, and a component of magnitude B±_ = Bq sin 6 
which appears to be rotating, with angular velocity ui relative to the grain. 
The rotating field B_l is perpendicular to u). 

It seems clear that the grain material will develop a magnetization M 
which will have a static component M|| = xoB|| parallel to cj, and a com- 
ponent M_L which rotates in body coordinates (but is stationary in inertial 
coordinates), with magnitude M± = \x{u))\B±, where xi^) is the (complex) 
magnetic susceptibility for a magnetic field rotating at frequency lu. In body 
coordinates, the rotating magnetization Mj, will "lag" behind the rotating 
field B_l - there will be an "in-phase" component xi('^)-B_l and an "out-of- 
phase" component X2{^)B±. 

The galactic magnetic field Bq will exert a torque on the grain magnetic 
moment VM. The torque is obviously _L to Bq. A bit of calculation shows 
that the out-of-phase component of the magnetization results in a torque of 
magnitude Vx2B±Bo = Vx2BQSm.9 which acts to leave the component of 
the angular momentum || to Bq unchanged (the torque, after all, must be 
_L Bo) but acts to reduce the component Jj_ = Jsin^ perpendicular to Bq: 



At low frequencies the imaginary component of the magnetic susceptibility 
varies linearly with frequency: X2 = Kuj. With w = J// = J\\/I cosO, noting 



-J^ = -Vx2Blsme 



(115) 



that J|| = const, we obtain 




tan^ 



(116) 



with solution 



tan 6* = tan(6'o)e-*/^°'^ 
I 



(117) 
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The magnetic susceptibilities of paramagnetic materials are quite well under- 
stood at the frequencies of interest (see, e.g., Morrish 1965), with 

If «1.25x 10-"sf^) . (119) 



With this value, and / w {2/5)pVa'^ , the Davis-Greenstein alignment time 



,2 



^3gcm-3 J \20KJ \ Bo J 

Thermally-rotating grains are disaligned by random collisions with gas atoms 
on a time scale Tdrag- The ratio of the alignment time to the disalignmcnt 



time is 



2 / „ X 1/2 



Tdrag ■ \iQcxar^ ) \2QK) \ Bq ) ^IOOK; VlQ-^cm^ 

(122) 

where we have assumed the hydrogen to be atomic. Thus we see that even for 
a = 10^^ cm grains and a relatively high value of the magnetic field strength, 
alignment is slow compared to disalignmcnt, and the ratio is linear in the 
grain radius a. The Davis-Greenstein alignment mechanism favors alignment 
of smaller grains, but this is a problem, since small grains are observed to 
NOT be aligned, while large grains are (Kim & Martin 1995). As we will see 
below, some important physics has been overlooked. 



7.5 Effect of Suprathermal Rotation 

Consider now a grain rotating suprathermally. The Davis-Greenstein align- 
ment time is independent of the rotation rate, since the torque is propor- 
tional to the rate of rotation. On the other hand, the time scale for random 
collisions to change the direction of the angular momentum is approximately 
{J/ Jth)^Tdriig, where Jth is the angular momentum expected for Brownian ro- 
tation. So if the grain rotation is suprathermal, with J/ Jth 3> 1, we can virtu- 
ally ignore the disaligning effects of random collisions, and Davis-Greenstein 
paramagnetic dissipation will bring the grain into alignment with Bq in a few 
million years. As we have seen in H6.2I suprathermal rotation appears to be 
expected. Could it for some reason be suppressed in small grains? 



7.6 Thermal Flipping vs. Suprathermal Rotation 

Why aren't small grains aligned? From eq. H97|) it certainly appears that H2 
formation should be able to drive grains as small as 10"^ cm into suprather- 
mal rotation, and yet there is strong evidence that grains smaller than 10~^ cm 
are not aligned (Kim & Martin 1995). 
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Lazarian & Draine (1999a) have proposed an explanation for this - the 
process of "thermal flipping" may prevent small grains from achieving suprather- 
mal rotation. 

The key issue is that the grain does NOT behave like a "rigid body" 
subject only to external torques - the grain has internal vibrational degrees of 
freedom. The grain rotational kinetic energy depends on the angle (3 between 
principal axis ai and J through eq. H104I) : 

E,,, = ^ + A^^—^s\n^P . (123) 

In the absence of external torques (and ignoring the small amount of angular 
momentum in the spin system), J is fixed, but the rotational kinetic energy 
can change if energy is exchanged with the vibrational degrees of freedom of 
the grain. Thus we expect the angle j3 to fluctuate. Suppose, for example, 
that we start at /3 = 0. Thermal fluctuations will cause (3 to "explore" the 
region < ^ < 7r/2. If 

J^iIl^lA<kT , (124) 

then thermal fluctuations could allow P to reach 7r/2: the maximum possi- 
ble rotational kinetic energy. At this point P can go either way, with equal 
probabilities, so the grain may "flip over" to the state tt/2 < (3 < tt. 

Now consider the effects of II2 formation torques. These torques are fixed 
in body coordinates. If the grain "flips" from /3 = to /? = tt, the torques 
change sign in inertial coordinates! As a result, if the grain flips frequently, 
the H2 torques can be averaged out to zero, and have no signiflcant effect on 
the grain dynamics, so that the grain rotation remains essentially thermal! 
Lazarian & Draine (1999a) refer to this as "thermal trapping". 

The question now comes down to estimating the rate of flipping, and 
how the rate depends on grain size. Lazarian & Draine (1999a) discussed 
the role which the Barnett effect (i.e., the electron spin system) can play 
in coupling the grain rotation to the vibrational modes, and concluded that 
flipping would be rapid enough to result in thermal trapping for grain sizes 
a^O.l/im. This seemed to successfully account for the absence of grain 
alignment for a ^ 0.1 /xm. 

However, further analysis of grain dissipational processes showed that 
the nuclear spin system could be important at the lower rotational rates 
appropriate to larger grains (Lazarian & Draine 1999b). It was found that 
the coupling of grain rotation to the nuclear spin system would lead to rapid 
thermal flipping even for grains as large as ~ 1 /im! 

Thermal flipping therefore now appears able to prevent II2 formation 
torques from being able to achieve superthermal rotation even for the a « 
0.2 fiiji grains which are observed to be aligned! How is grain alignment to be 
achieved? 
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7.7 Effect of Radiative Torques 

So long as the grain surface is not modified, the H2 formation torques are fixed 
in body coordinates, and thermal flipping will cause these torques to average 
to zero if the grain spends equal amounts of time in the two "flip states" 
(i.e., /3 < 7r/2 and /3 > 7r/2). However, the grain is also subject to other 
systematic torqiies, the most important of which is the radiative torque due 
to anisotropic starlight. This torque is not fixed in body coordinates, so that 
even if the grain spends equal amounts of time in the two flip states, there is a 
secular change in angular momentum due to the radiative torques. Radiative 
torques are unimportant for small a ^ 0.05 /um grains, but very important for 
a ^ 0.1 /zm grains. It therefore appears that "thermal trapping" will prevent 
small (a ^ 0.05 /im) grains from achieving suprathermal rotation, but that 
larger (a ^ 0.1 /um) grains can achieve suprathermal rotation, following which 
they can be gradually aligned by the Davis-Greenstein mechanism. 

But life turns out to be more complicated. Draine & Weingartner (1997) 
found that radiative torques act not only to spin up grains, but also to change 
their alignment. The characteristic time scale for radiative torques to change 
the grain alignment turns out to be just the drag time Tdrag, which is shorter 
than the Davis-Greenstein alignment time, so that the radiative torques ap- 
pear able to dominate over the systematic effects of the Davis-Greenstein 
alignment torque. 

Study of three different grain shapes and different angles between the 
magnetic field direction and the starlight anisotropy direction showed that 
under some conditions the starlight torques brought about alignment (the 
aligned state was an "attractor"), while under other conditions the grain did 
not settle down in any stationary state. There appeared to be an overall 
tendency for grain alignment to take place, but the sample of 3 grain shapes 
was insufficient to reach general conclusions. 

Further study of this problem is underway, and it is hoped that in a year 
or two we will have definite conclusions on the role of radiative torques in 
the grain alignment problem. But for the moment the situation appears very 
promising - we seem to finally understand the absence of alignment of small 
grains, and it appears that the physics we now understand will explain the 
alignment of the larger ones, with radiative torques due to starlight playing 
a critical role. 

7.8 Alignment and Disalignment in Different Regions 

There is evidence that the degree of grain alignment varies from one region to 
another. Observations of polarization of starlight indicate that the a ^ 0. 1 fim 
grains in diffuse clouds are substantially aligned, but JHK observations of 
quiescent dark clouds (Goodman et al. 1995) indic;ate that the dust grains in 
the dark interiors are not aligned. How can we understand this? 
There are several possible plausible answers: 
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Fig. 37. 60 /im linear polarization toward M17. Flux contours are at 20%, 30%, 
...,90% of the peak flux. From Dotson et al (2000). 



• The inner dark regions have insufficient starhght for the radiative torques 
to achieve grain alignment. 

• The inner dark regions are molecular, so that H2 formation torques can- 
not drive suprathermal rotation. 

• Photoelectric emission is suppressed in the inner dark regions. 

It is also possible, of course, that either the magnetic field structure is com- 
pletely tangled in the interior, or the grains in the central regions are for 
some reason more spherical, but these do not seem plausible. 

It is also interesting to note that in some dark clouds the grains are rea- 
sonably well aligned: for example, the M17 molecular cloud fFig. I37|) shows 
linear polarization (averaged over the 22" beam) as large as 5%. Note, though, 
that although this cloud is "dark" (i.e., we cannot see into it at optical wave- 
lengths) it is full of stars - the infrared spectrum indicates a grain tempera- 
ture T sa 45 K, and the grains must be heated by starlight: the mean starlight 
intensity must be ^ (45/20)^ w 10^ times greater than the ISRF in order to 
heat the grains to this temperature! It therefore appears likely that starlight 
torques could drive the grains in this cloud to suprathermal rotation rates. 
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Finally, a recent study by Padoan et al. (2001) concluded that submm 

polarization maps of quiescent cores were consistent with a model where the 
grains in the outer parts of the cloud with Ay ^ 3 mag were aligned, while 
grains in regions with Ay > 3 mag were randomly oriented. 

All of the above appear consistent with the idea that starlight torques 
are an essential part of the grain alignment process. Grains deep in quiescent 
dark clouds are unable to align because starlight torques are weak or absent, 
whereas starlight torques are present are do produce alignment within star- 
forming clouds. 

7.9 Summary of the Grain Alignment Story 

The grain alignment story has been long and complicated, with a number of 

false turns. The discussion above has skipped some of the important physics 
involving coupling of the grain rotation to the electron and nuclear spin 
systems (see Lazarian & Draine 1999a,b) as well as the important question 
of the grain dynamics during "crossovers" , when reversals of the systematic 
torques on the grain cause the grain angular momentum to become small 
- during this time the grain is very vulnerable to disalignment (see, e.g., 
Lazarian & Draine 1997). 

To summarize: the important elements of the grain alignment story now 
appear to be as follows: 

• Grains are subject to 2 important systematic torques: 

— H2 formation torques, 

— radiative torques due to starlight. 

• For grains with rotational kinetic energies less than a few times kTd, 
thermal fluctuations allow the grain to sample different orientations of 
the principal axes relative to the instantaneous angular momentum J, 
leading to "thermal flipping" of the grain between the two "flip states" . 

• If the thermal flipping is sufflciently rapid, the grain will spend 50% of 
the time in each flip state, in which case the H2 formation torque, which 
is flxed in body coordinates, will integrate to zero in inertial coordinates: 
the H2 formation torque will be unable to drive the grain to suprathermal 
rotation. 

• Starlight torques are not fixed in body coordinates, and hence do not 
integrate to zero even if the grain spends equal amounts of time in the 
two flip states. 

• For grains with a ^ 0.1 fxm in regions with sufficiently strong starlight, 
the systematic torque due to starlight will be able to drive the grain to 
suprathermal rotation. 

• Once a grain achieves suprathermal rotation, thermal fiipping will cease, 
the grain will settle into a state where the principal axis ai is cithcir 
parallel or antiparallel to J, and all systematic torques will operate si- 
multaneously. 
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• Grains in suprathermal rotation are essentially immune from disalign- 
ment by random collisions. Alignment will result from the combined ef- 
fects of Davis- Greenstein alignment and starlight torques. 

• Under some circumstances, starlight torques do not result in strong align- 
ment. However, numerical experiments suggest that, averaged over irreg- 
ular grain shapes, starlight torques lead to overall alignment. 

• The fact that starlight torques are important only for grains with a ^ 0.1 yum 
appears to account for the observation that a ^ 0.05 /im grains in diffuse 
clouds are not aligned, while a ^ 0.1 /xm grains appear to be well-aligned. 

8 Evolution of the Grain Population 

The interstellar grain population in a galaxy is the result of a rich and com- 
plicated mix of processes. 

We would like to know the composition of interstellar dust, how and where 
it forms, and the extent to which individual grains are either homogeneous or 
conglomerate. We would like to know what the size distribution of the differ- 
ent types of grains is, how and where the size distribution is fashioned, and 
the manner in which the size distribution varies from one region to another. 
If we understood the balance of processes acting in our part of the Milky 
Way, perhaps we could deduce how grains might differ in the inner or outer 
parts of the Milky Way, and in the LMC, SMC, and other galaxies. 

Unfortunately, the overall grain evolution problem seems to me to be 
beyond our grasp at this time I don't think we understand the parts well 
enough to justify trying to put them together as a whole. At this time, we 
have to try to clarify the individual elements of the story. 

8.1 Dust Formation in Stellar Outflows 

There is abundant evidence that dust grains form in some stellar outflows: we 
see infrared emission from dust around red giants, around carbon stars, and 
in planetary nebulae. In all these cases the dust condensed out of material 
which at an earlier time was completely gaseous. Observations show that 
different types of stellar outflows form different types of dust. 

The dust around red giants with 0/C > 1 shows a 10 /xm emission feature 
which is due to the Si-0 stretching mode in amorphous silicate material; 
around some OH/IR stars the 10 /xm feature appears in absorption (e.g., 
Demyk et al. 2000), along with emission features revealing the presence of 
crystalline silicates. 

The dust around carbon stars, on the other hand, does not show the 10 /xm 
feature, which is consistent with the fact that we do not expect silicates to 
form in an atmosphere with 0/C < 1. Some carbon stars (e.g. IRC+10216) 
instead show a 11.3 /itm emission feature which is attributed to SiC (e.g., 
Blanco et al. 1998). 
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It should also be mentioned that we find presolar grains in meteorites 
whose isotopic composition clearly indicates formation in outflows from evolved 
stars (Hoppe & Zinner 2000); these include nanodiamonds, SiC, graphite, and 
AI2O3. Some of these grains may have formed in supernova ejecta. 

So there is no question that stellar outflows contribute dust to the inter- 
stellar medium. However, this does not necessarily imply that the bulk of the 
dust in the interstellar medium was formed in stellar outflows. 

The physics and chemistry of dust formation in stellar outflows is com- 
plex: the material is generally far from thermodynamic equilibrium, and the 
outflows themselves are often (perhaps always?) hydrodynamically complex 
- neither steady nor spherically symmetric. We can not yet reliably describe 
the details of the grain formation process. 

8.2 Grain Grovifth in the ISM 

lOl l Accretion in the ISM 

Dust grains in the interstellar medium can grow, since atoms from the 
gas can stick to them. We can calculate the lifetime against accretion for an 
atom or ion in the gas phase: 

= A7'/'s. ( ^) / dana'^Dia) , (125) 
V™H/ J da 

where Ai is the mass number of the ion, Si is the sticking coefficient, and nd{a) 
is the number density of grains with radii less than a. The enhancement factor 
D{a) gives the collision rate relative to what it would be in the absence of 
electrostatic effects: 

D{a)^Y.f^^d,a)B{Zd,a) , (126) 

where 

{exp {-ZdZ,e^/akT) for ZgZ, > 
(1-^)^ forZ,Z.<G (^27) 
forZ, = 

Fig. 1381 shows the collisional enhancement factor D{a) calculated by Wein- 
gartner & Draine (1999) for carbonaceous and silicate grains in various envi- 
ronments. Enhancement factors D ^ ^ 10^ are found for a ^ 10^^ cm grains 
in the CNM. As a result, an ion (e.g.. Sill, Till, Fell) would have a lifetime 
against accretion of ~ 1 x 10^ yr in the CNM if the grain size distribution 
includes the ultrasmall grain population of Fig. EH Therefore we expect that 
those metals which can "stick" to the small grains will be strongly depleted 
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Fig. 38. Collision enhancement factor for positively charged ion (see text). From 
Weingartner & Draine (1999). 

from the gas phase; presumably those elements which are not strongly de- 
pleted (e.g., Na, K, S) for some reason do not remain on the ultrasmall grains 
after they collide. 

101 2 Mantle Formation 

As discussed in ill. 81 grains in dark molecular clouds are apparently coated 
with ice "mantles". The overall mantle composition is dominated by H2O ice. 
It is not known whether the H2O molecules arrived from the gas phase and 
then froze onto the grain, or whether the H2O formed by surface chemistry. 
In most clouds the gas phase II2O abundance is very low, and the observed 
mantles would not be able to form by accretion, suggesting that O and H 
atoms arriving at the grain react there to form H2O (Jones & Williams 1984). 
However, it has recently been proposed that the II2O could be produced in 
the gas phase during occasional passage of a shock wave, and then quickly 
accreted onto the grain surfaces (Bergin et al. 1999). Other molecules (e.g., 
CO) may be mixed with the H2O, or in some cases condensed separately, 
as the absorption profile in some cases requires that the CO not be mixed 
with polar molecules. This could result if the H2O was deposited while the 
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grain was too warm to retain CO, with the CO later freezing out when the 
conditions changed and the grain became colder. 

101 3 Coagulation 

The grain population in normal diffuse clouds has a very large population 
of very small grains, required to explain the rapid rise in extinction at short 
wavelengths. In dense regions, however, the extinction at short wavelengths 
appears to be reduced. This can only happen if some of the small grains are 
either destroyed altogether, or simply coagulated onto big grains. The latter 
seems more likely. 

What is the time scale for such coagulation? This depends on the velocity 
differences between different grains. Grain-grain velocities can be the result 
of differential drift velocities (see ^5.5|l . hydrodynamic turbulence (see, e.g., 
Draine 1985), or MHD turbulence (Lazarian & Yan 2002). 

Grains larger than ~ 10~^ cm contribute a geometric cross section per H 
E Ri 0.5 X 10^^^ cm^. If the characteristic grain-grain velocity difference is 
Avci, then the grain-grain coagulation time would be 

r ^.1-1 „ /30cm^3\ /0.1kms-i\ 

tcoag = [niiSAvd] ^ = 2x 10^ yvi——j , (128) 

so that grain coagulation will not be very important in diffuse clouds, but 
could proceed in dense regions. 

8.3 Photolysis 

When grains with ice mantles are returned to diffuse regions and exposed to 
ultraviolet radiation, it appears that most of the material is driven off, but 
the ultraviolet radiation photolyzes the ice to produce an organic refractory 
residue (e.g., Greenberg et al. 2000). This residue may be responsible for the 
interstellar 3.4 fim "aliphatic C-H" absorption band. 

8.4 Photodesorption 

The process of photodesorption is potentially very important, but is unfortu- 
nately not well understood because of limited laboratory studies. Photodes- 
orption is a quantum process whereby absorption of a single photon can lead 
to ejection of an atom or molecule from a solid. For example, an absorbed 
photon might cause an atom or molecule on a substrate to be excited to a 
repulsive electronic state which would then be ejected from the substrate. 

Photodesorption was discussed by Draine & Salpeter (1979b), who argued 
that a molecule in a surface monolayer might have a photodesorption cross 
section as large as ~ 10^^^ cm^. 

If an atom or molecule adsorbed on a grain surface has a cross section 
apd = 10~^*cr_i8 cm^ for S < hv < 13.6 eV photons, the photodesorption rate 
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for an adsorbed atom or molecule would be 6 x lO~^^0-i8 in the diffuse 
ISM, since the density of 8 — 13.6 cV photons is ^ 2.0 x 10~"^cm~'^ in the 
MMP ISRF. Laboratory studies show that H2O ice has a photodesorption 
cross section apd^Sx 10~^^ cm^ at A = 1215A (Westley et al. 1995a, 1995b), 
so that it seems likely that at least some species will have (T_i8 ^ 1, with pho- 
todesorption rate ^ 6 x 10~^^ s~^. Note that this estimated photodesorption 
rate is comparable to photodissociation rates for diatomic molecules in the 
interstellar radiation field. ^ 

O is the most abundant species after H and He. The probability per unit 
time that an O atom will arrive at a single surface site (of area ^ 10~^^ cm^) 
is ^ 2 X 10~^^s^^ in the CNM. It is therefore apparent that any species 
with (7i8 ~ .01 will be photodesorbed before accretion of anything chemically 
interesting (other than H) can take place on top of it. 

We conclude that photodesorption most likely plays a dominant role in 
determining what kind of accretional growth takes place on a grain. It is no 
surprise that the noble gases (He, Ne, Ar) do not deplete, but the fact that 
some chemically-reactive elements (e.g., Na, K, S) appear to undergo minimal 
depletion, at least in the diffuse ISM, may be due to photodesorption.^ Lab- 
oratory studies of photodesorption from carbonaceous or silicate substrates 
would be of great value to clarify why some elements deplete and others do 
not in the diffuse ISM. 

8.5 Grain Destruction in Shock Waves 

From time to time a shock wave will pass through an clement of interstellar 
gas, compressing, heating, and accelerating the gas (Draine & McKee 1993), 
and creating conditions under which grain destruction can occur (Draine 
1995). Shock waves can be driven into cold cloud material ahead of an ion- 
ization/dissociation front (Bertoldi & Draine 1996), as the result of energetic 
outflows from newly- formed stars, or as the result of the explosion of a nearby 
supernova. 

It is easiest to study the shock in the "shock frame" where the shock front 
is stationary: in the upstream direction the grains and gas move together 
toward the shock with a speed equal to the shock speed Vs . The fluid is then 
decelerated at the shock front, with the postshock gas velocity Vg ^ for a 
strong shock. The postshock gas is initially hot. If the shock is strong, and we 
can neglect the effects of a possible magnetic precursor, the gas temperature 
will rise to 

^ Roberge et al. (1991) estimate photodissociation rates 8.6 x 10 s ^ and 3.5 x 
10"^° s"^ for CH and OH, respectively. 

^ The only other depletion-suppressing process would be chemisputtering by H 
atoms: an impinging H atom might react with the adsorbed atom, and the reac- 
tion energy might eject the resulting hydride. 
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where /i is the mass per free particle {ii/mn k, 1.4/1.1 = 1.27 for cold atomic 
clouds, and ii/mn = 1.4/0.6 = 2.33 for molecular gas, where we assume 
He/H=0.1). 

The grains move across the shock front into the postshock fluid, where 
they initially have velocity Vs — Vg (3/4)ws relative to the gas. The grains 
now begin to decelerate by gas drag, but are also acted on by the magnetic 
field B and electric field E — — x B/c (evaluated in the shock frame). As 
a result, the grain population develops grain-grain velocity differences which 
can be of order the shock speed Vs , or even larger if conditions are favorable 
for "betatron" acceleration of the dust as the density and magnetic field 
increase in cooling postshock gas (Spitzer 1976; Draine & Salpeter 1979b). 

18.51 1 Grain-Grain Collisions 

In the absence of betatron acceleration, a grain will slow down by a factor 
1/e when it has collided with its own mass of gas. If all grains were identical, 
then the chance that a grain would run into another grain before appreciable 
slowing-down would be of order the dust-to-gas mass ratio of ~ 0.01, and we 
can therefore expect that a typical grain (one representative of the bulk of 
the grain mass) will have a probability of ~ 0.01 of collision with another 
comparable grain with a grain-grain velocity difference ~ Vs . 

The physics of grain-grain collisions has been discussed by Tielens et al. 
(1994). For refractory grains one expects virtually complete destruction if the 
energy in the center-of-mass system is sufficient to vaporize both grains. For 
identical grains, mean atomic mass 20 amu, and binding energy per atom 
of 5.7 eV (these numbers are appropriate to MgFeSi04 silicate) complete 
vaporization could occur in a head-on collision between identical grains with 
a velocity difference of 15 kms~^, so we may expect vaporization of perhaps of 
1% of the grains in a hydrodynamic shock with shock speed Vs ^ 20km s^^. 
At lower grain-grain collision speeds, shattering may occur. Thus we should 
expect that several percent of the large dust grains will be shattered in shocks 
with Us XL 10km s^^. Grain-grain collisions may be responsible for maintaining 
the population of small grains in the interstellar medium. 

1^31 2 Sputtering 

In a fast shock, the grain material can be eroded by the process of "sput- 
tering", where individual atoms or ions from the gas collide with the grain, 
occasionally ejecting one of the grain atoms. In interstellar gas, the sputtering 
rate depends mainly on the sputtering "yields" Y{E) for H and He projec- 
tiles incident on silicate or carbonaceous target materials. In ionized gas the 
sputtering rates will be affected by the grain charge. An expression for the 
rate of sputtering of a moving grain, and estimates for sputtering yields Y , 
are given by Draine & Salpeter (1979a). 

In the initial postshock region where the gas is hot, the sputtering is due 
to a combination of thermal energy and the motion of the grain through the 
gas with velocity ~ Svg/A. If the gas cools before the grain slows down, then 




Fig. 39. Sputtering yields for H and He on carbon estimated by Barlow (1978), 
Draine & Salpeter (1979a), and Tielens et al. (1994), and compared with laboratory 
data. 



subsequent sputtering is entirely due to the rapid motion of the grain through 
the gas, with He atoms having kinetic energies 4 times larger than H atoms. 
The sputtering continues until the grain slows to the point where the kinetic 
energy of the impinging He atoms is too low for appreciable sputtering. 

Detailed models of grain destruction in shock waves are in general agree- 
ment that shock speeds Vs ^ 100km s^^ result in substantial grain destruc- 
tion, with essentially complete grain destruction for Vs > 200km s~^ in the 
case of a radiative shock - see Fig. 1411 (Draine & Salpeter 1979b; Jones et al. 
1994). 

These estimates for grain destruction sppear to be in agreement with the 
observation that high-velocity gas (which has been shocked) tends to have 
enhanced gas phase abundances of elements, such as Si, which are normally 
heavily depleted in interstellar gas, an effect first noted for Ca by Routly & 
Spitzer (1952). Fig. shows this effect for Si. 

8.6 Thermal Sublimation 

Ice mantles on grains can be removed if the grains become too warm. What 
is the critical grain temperature for ice mantle removal? 




Fig. 40. Sputtering yields for H and He on silicate material estimated by Barlow 
(1978), Drains & Salpeter (1979a), and Tielens et al. (1994). 




V. ( km s ' ) ( km s ' ) 



Fig. 41. Fractional grain destruction calculated for graphite and silicate grains in 
interstellar shock waves in low density gas. Figure from Jones et al. (1994). 



The probability per unit time for a surface molecule to desorb is approx- 
imately 

at 

where f/o is the binding energy of the molecule, and i/o ~ lO^'' s^^ is a char- 
acteristic vibration frequency. If there are iYmono monolayers of ice on the 
grain and we want it to survive a time At, the grain temperature should not 
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Fig. 42. Observed gas-phase Si abundances: evidence for refractory grain destruc- 
tion in high-velocity shocks. Solid, dashed, and broken lines are theoretical models 
discussed by Cowie (1978). Figure from Cowie (1978). 



exceed 

Uo/k 

'^"•^ - ln(z.oZ\t/iV„.ono) ^ ^ 

Some heating phenomena are brief (e.g., a stellar flare or a supernova) but 
we are usually interested in longer time scales. For example, a massive star 
in or near the cloud might have a lifetime At h x 10^ yr. We can write 

- - ^^/^^^ (132) 



1 + 0.018 ln[(Z\t/106 yr)(10ViVmono)] 

Thus, a H2O ice mantle (C/q = 0.5 eV) 10"^ monolayers thick will survive 
10^ yr provided T < 108 K. Therefore we should not be surprised to find 
H2O ice in clouds where the grain temperature is 80 K (as in the molecular 
cloud OMC-1), but it would be a great surprise to find H2O ice if the grain 
temperature were, say, 120K. 



8.7 Gamma Ray Bursts 

Gamma ray bursts may destroy dust grains out to a substantial distance. 
While this is probably uninteresting in terms of evolution of the overall grain 
population in galaxies, it may nevertheless be of observational interests be- 
cause it affects the observable properties of gamma-ray bursts. This is of 
considerable interest because there is evidence that gamma-ray bursts are 
associated with star-forming regions in galaxies, and yet the optical after- 
glows from a number of bursts show little evidence for reddening. 

Grain destruction near GRBs has been discussed by Waxman & Draine 
(2000), Fruchter, Krolik & Rhoads (2001), and Draine & Hao (2002). There 
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are two separate mechanisms which could be important. Thermal sublimation 
(as discussed above, but now with larger binding energies Uo/k ~ 7 x 10'' K 
and shorter times At « 10 s for the optical transient, so the grain temperature 
can rise to ~ 3000 K. Draine & Hao find that for a reasonable (but very 
uncertain) estimate of the GRB peak optical luminosity, dust grains can be 
destroyed out to distances of ^--^5 pc. 

A more exotic grain destruction mechanism may operate; if effective, it 
could destroy grains out to greater differences. This is the mechanism of 
"Coulomb explosion" (Waxman & Draine 2000). The grains will be quite 
highly charged due to the intense X-ray irradiation from the GRB. If the 
grains are sufficiently weak, the electrostatic stresses on the grain might con- 
ceivably tear it apart. Fruchter et al give arguments in favor of this outcome, 
in which case GRBs could destroy grains out to quite large distances. 

8.8 Time Scales 

When all grain destruction processes are considered, it is found that sput- 
tering in shock waves dominates the overall destruction of refractory grains. 
The frequency with which grains are overtaken by such shocks depends on 
the frequency of the supernovae which drive these blastwaves, and on the 
structure of the interstellar medium, since material in dense regions is "shel- 
tered" from the effects of a passing blastwave. The problem is complex, and 
therefore there are many places where different assumptions could alter the 
quantitative conclusion, but Draine & Salpeter (1979b) estimated a lifetime 
^ 1.3 X 10^ yr for silicate material in a diffuse cloud, while Jones et al. (1994) 
find a lifetime 2.2 x 10* yr for silicate grains. 

Draine & Salpeter (1979b) stress that this relatively short lifetime was 
impossible to reconcile with the idea that "interstellar grains are Stardust" : 
nearly all of the Si in the interstellar medium is in solid form (e.g., see Fig.|51 
showing that only 2.5% of the Si is in the gas phase in the cloud toward COph). 
The "turnover" time of the interstellar medium is equal to the ratio of (mass 
in gas)/ (rate of mass going into stars), or ~ 5xlO®Mo/5M0yr~'^ « IxlO^yr. 
Suppose that all the Si leaving stars entered the ISM in solid form, but there 
was no grain growth in the ISM. A grain destruction time ^ 2 x 10* yr - 
only ^-^20% of the turnover time - would imply that ^ 80% of the Si would 
be found in the gas phase. However, observations of diffuse clouds typically 
show < 10% of the Si in the gas gas phase (see, e.g.. Fig. 1^1). 

The ineluctable conclusion is that most of the Si atoms in grains must 
have joined those grains in the ISM through some process of grain growth. 
The kinetics of accretion implies that only in denser regions can accretion 
rates be large enough to remove Si from the gas fast enough to maintain the 
low observed gas-phase abundances. This in turn requires that there must 
be fairly rapid exchange of mass between the dense and less dense phases of 
the ISM. A more quantitative discussion of the required mass exchange rates 
may be found in Draine (1990) and Weingartner & Draine (1999). 
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9 Effects of Dust on Interstellar Gas 

9.1 Photoelectric Heating 

Photoelectrons emitted from dust grains depart with nonzero kinetic energy. 
This is often the dominant mechanism for heating intersteUar gas, whether 
in diffuse clouds or in photodissociation regions. 

Photoelectric heating rates have been reestimated recently by Bakes & 
Tielens (1994) and Weingartner & Draine (2001c). For starlight with the 
spectrum of the ISRF, Weingartner & Draine find a heating rate per H nu- 
cleon r/n-a « 6 x 10~^^ erg s~^ for conditions characteristic of cool H I clouds. 
This heating rate is much larger than other heating mechanisms, such as ion- 
ization by cosmic rays or X-rays. 

9.2 H2 Formation and Other Chemistry 

The H2 molecule is central to interstellar chemistry. The rovibrational lines 
can be an important coolant for gas at temperatures 100 ^ T ^ 2000 K. 
The principal gas-phase formation processes are via H~ 

ll + e-->ll-+hiy ; + H ^ H2 + (133) 

or 3-body collisions: 

H + H + X^Ha + X + Zii; (134) 

where the third body X could be H, He, or another H2. The H~ formation 
pathway H133|l is slow because of the generally low abundance of H~ in diffuse 
gas, and the 3 body process (jl34|l requires very high densities to become 
important. In the interstellar medium of the Milky Way, the dust abundance 
is such that H2 formation is overwhelmingly due to grain surface catalysis. 
The grain catalysis process was first proposed by Gould & Salpeter (1963), 
and the basic picture has not changed: 

1. An H atom arrives at a grain surface and "sticks". 

2. The H atom diffuses over the grain surface until it becomes trapped by 
either chemisorption or physisorption 

3. A second H atom arrives at the surface, sticks, diffuses over the surface, 
and encounters the trapped H atom. 

4. The two H atoms react to form H2, releasing ~ 4.5 eV of energy, and 
ejecting the H2 from the grain surface into the gas. 

Part of the 4.5eV goes into overcoming the binding of the two H atoms to 
the surface, part goes into vibrational excitation of the grain lattice ("heat"), 
part appears as kinetic energy of the H2 , and part appears in vibrational and 
rotational energy of the newly-formed H2 molecule. 
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In the above picture, the H2 formation rate is 

rfn(H2) 



dt 

1/2 



Rniin{B.) (135) 



^=2 



fda^'^.aHa) (136) 
2 \TTmiiJ J nil da ' ^ ' 

7.3 X 10-^ cm^s-^(-Z^) ,._2f„^lL-i (137) 



^lOOKy 10-21 cm2H-i 
A = y da(l/nYi){dnd/da)iTa'^ (138) 

{l) = \j dail/nH){dnd/da)wa^^{a) (139) 

where 7(a) is the fraction of H atoms coUiding with a grain of radius a which 
are converted to H2; 7 < 1 since not all impinging H atoms will "stick", and 
some of those which stick may later be removed from the surface by some 
process other than H2 formation. Gould & Salpeter argued for 7 of order 
unity under interstellar conditions. 

Dust models which reproduce the observed extinction of starlight have 

A= [ da—^na^ > lO'^^ cm^ (140) 

J «H da 

For example, the dust model of Weingartner & Draine (2001a) has A = 
6.7 X 10-21 cm2/H. The H2 formation "rate coefficient" R can be determined 
from ultraviolet observations of H2 in diffuse regions; Jura (1975) found R Ri 
3 X 10-^^cm^s-^, which implies (7) « 0.06. It appears, then, that the very 
small grains which dominate the grain surface area (in particular, the PAHs) 
must not be efficient catalysts for H2 formation. 

How might H2 formation on the smallest grains be suppressed? The time 
scale between photon absorption events is 

Tabs ~ 10S(10A/a)3 , (141) 
while the time scale between arrival of H atoms is 

.H«8xl0«s(l^)'(^l fi^'^'^ 



(142) 



so it is possible that an H atom physisorbed on a very small grain may be 
removed when the grain is heated by a photon absorption event, so that 
only rarely would an arriving H atom find another H atom with which to 
recombine. The efficacy of this process depends on the binding energy of the 
H atom. 
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9.3 Grains as Sources of Complex Molecules 

Observational studies of star-forming regions sometimes find molecular abun- 
dances in the gas phase which arc difhcult to understand in the context of 
pure gas-phase chemistry. An extreme example of this is the observation 
of D2CO/H2CO=0.003 toward Orion (Turner 1990), or D2CO/H2CO=0.03- 
0.16 toward the low-mass protostar IRAS 16293-2422 (CeccarelH et al 2001). 
Such extreme deuteration seems impossible to envision in the gas phase, but 
could occur by chemistry on grain mantles during a precollapse phase, with 
the molecules put into the gas phase when the dust is warmed by energy 
from the protostar. There seems little doubt that at least some regions have 
molecular abundances which are heavily influenced by grain surface chem- 
istry. 

9.4 Ion Recombination on Dust Grains 

In diffuse clouds, the ion/neutral fraction elements with ionization potentials 
< 13.6 eV is indicative of the relative rates for photoionization of starlight vs. 
recombination of ions with electrons. It is often assumed that the only channel 
for neutralization of metal ions X+ is radiative recombination, X+ + e~ — > 
X + hv. However, collisions with neutral or negatively charged grains can 
be more important than radiative recombination for neutralization of metal 
ions for conditions typical of the diffuse interstellar medium (Weingartner & 
Draine 2001d). Weingartner & Draine provide estimates for the effective "rate 
coefficient" for ion neutralization via collisions with dust grains. Recombina- 
tion on dust grains is particularly effective for protons in H I regions, and 
therefore is involved in regulating the free electron density. 

9.5 Coupling Neutral Gas to Magnetic Fields 

In gas of low fractional ionization, a significant fraction of the "free" charge 
present may be located on dust grains. The fraction of the free charge residing 
on the grains depends very much on the numbers of very small dust grains, 
but appears likely to become significant in regions of low fractional ionization 
ne/nn ^ 10"^ (Draine & Sutin 1987). 

Coupling of the magnetic field to the gas is due to the force exerted on 
moving charges; in regions of low fractional ionization, the neutral gas will be 
unaffected by the magnetic field unless the charged species are coUisionally 
coupled to the neutral atoms and molecules. Because dust grains have rela- 
tively large physical cross sections, the charge trapped on dust grains can be 
important at coupling the magnetic field to the neutral gas in MHD shocks in 
gas with fractional ionization rie/nn ~ 10"^ (Draine 1980; Draine, Roberge 
& Dalgarno 1983). 
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9.6 Dust Grains as Magnetometers? 

Magnetic fields are dynamically important in the ISM, but difficult to observe 
remotely. Nature has been kind enough to strew microscopic magnetometers 
- dust grains - throughout the interstellar medium. Unfortunately we have 
not yet figured out how these magnetometers work, so that even if we could 
measure their degree of alignment, we would not be able to determine the 
magnetic field strength. But we can use dust grains as "compasses" to indicate 
the direction of the magnetic field. 

As discussed above, dust grains have their angular momenta J systemat- 
ically aligned with the magnetic field Bq, at least in regions where sufficient 
starlight is present. Here we review the main features of this alignment: 

• It is virtually certain that the sense of alignment is to have J || Bq (to 
within a sign - parallel and antiparallel are equally favored) . 

• Suppose the grain has 3 principal axes ai, 3.2, and a^, with moments of 
inertia Ii > I2 ^ la- If the grain has rotational kinetic energy i?i-ot > 
(3/2)fcTgi., then the grain will tend to have its shortest axis ai || J so as 
to minimize its rotational kinetic energy at constant J. In regions where 
the grains are cooler than the gas (the most common situation), this is 
expected to be true for both large grains (which rotate suprathermally) 
and small grains (which rotate thermally). 

• Thus we expect the grains to have their "long axes" (a2, £13) _L Bq. 

• Starlight extinction will be greatest for starlight with E || the long axes, 
and therefore we expect the transmitted starlight to have E || Bp. 

• Far-infrared and submm emission will be polarized |j to the grains' long 
axes, and hence we expect the FIR/submm emission to be polarized with 
E_L Bq. 

Therefore we can obtain a map of the magnetic field by either measuring 
the polarization of starlight for many stars, or by mapping the polarized 
far-infrared or submm emission from the cloud. CCDs and infrared detectors 
make possible extensive measurements of starlight polarization, even in dense 
regions, and FIR/submm arrays (e.g., SCUBA on the JCMT) make possible 
measurements of polarized submm emission. 

Grains are likely to have a fairly high polarizing efficiency in the FIR/submm. 
Define 

Cpoi = i [C(E II a2) + C(E II ag)] - C(E || ai) (143) 

Cran = I [C(E II ai) + C(E II a2) + C(E || a^)] (144) 

If the grains were spherical, there would of course be no polarization in ex- 
tinction or emission. Fig. 1431 shows how the polarizing efficacy depends on 
the axial ratio of the dust grains. The grains are assumed to be either oblate 
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Fig. 43. Polarization efficiency Cpoi/Cran for spheroidal dust grains (oblate and 
prolate) as a function of the ax;ial ratio. From Padoan et al. (2001). 

or prolate spheroids, and the dielectric function of the grains is taken to be 

that of "astronomical silicate" at 350 /im. 

For perfect spinning alignment, the maximum polarization is obtained for 
a sightline _L Bq: 



Suppose the grain shape has (intermediate+long)/(2 x short) = 1.4 (e.g., ax- 
ial ratios 1:1.2:1.6). Then Cpoi/Cran ~ 0.5 for astronomical silicate. Perfect 

alignment would give polarization in emission P = 0.5/(1 — .08) = 0.55. 
However, the largest observed polarizations are ~ 10%, implying either 

• imperfect alignment 

• less elongation 

• disorder in the magnetic field on short length scales. 

Unfortunately, it will be difficult to determine the relative importance of 
these three effects. For perfectly aligned grains in a uniform magnetic field 
perpendicular to the line-of-sight, P = 0.1 would be produced by grains with 
axial ratio of only ~ 1.05 - very minimal elongation. 



max — 




(145) 
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10 Concluding Remarks 

There are many unanswered questions in the astrophysics of interstehar dust 
grains in cold clouds. The questions are interesting in their own right, but 
one is also driven to answer them by the need to understand dust grains well 
enough to use them as diagnostics of interstellar conditions, and to under- 
stand the effects which dust grains have on the dynamics and evolution of 
interstellar gas. With the advent of powerful infrared and submillimcter ob- 
serving facilities - including SIRTF, SOFIA, ALMA, Herschel, and Planck 
- we can anticipate that there will be progress on the astrophysics of dust, 
driven by the challenge of understanding the new data. An exciting decade 
lies before us! 

I thank Robert Lupton for making available the SM software package. 
This work was supported in part by NSF grant AST 9988126. 
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